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FOREWORD 


This  report  is  one  of  a  series  of  related  papers  covering  various 
aspects  of  a  broad  program  to  investigate  the  flow-field  variables 
associated  with  hypersonic- velocity  projectiles  in  free  flight  under 
controlled  environmental  condtions.  This  research  is  being  conducted 
in  the  Flight  Physics  Range  of  General  Motors  Defense  Research 
Laboratories;  the  work  is  supported  by  the  Advanced  Research  Pro¬ 
jects  Agency  under  Contract  No.  DA-04-495-ORD-3567  (Z),  and  the 
contract  is  monitored  by  the  Army  Missile  Command.  It  is  intended 
that  the  series  of  reports  generated  under  this  program,  when  com¬ 
pleted,  shall  form  a  background  of  knowledge  of  the  phenomena  involved 
in  the  basic  study  and  thus  aid  in  a  better  understanding  of  the  data 
obtained  in  the  investigation. 

The  reports  submitted  under  this  program  fall  into  the  following 
general  categories: 

1.  Technical  discussions  of  the  theoretical  and  experimental 
aspects  of  selected  areas  of  interest, 

2.  Data  reports  and  compilations  of  experimental  data  derived 
in  the  investigation,  and 

3.  Periodic  progress  reports. 

The  present  report  is  one  of  those  in  the  first  category. 
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ABSTRACT 


The  usual  manifestation  of  the  phenomenon  known  as  the  radar 
absorption  effect  is  the  drastic  reduction  of  the  radar  cross  section 
of  a  body  when  it  is  surrounded  by  a  plasma  sheath  having  certain 
spatial  configuration  and  electromagnetic  properties.  This  effect  has 
been  observed  during  the  atmospheric  reentry  of  hypersonic  test 
vehicles  and  in  laboratory  studies  of  models  in  hypersonic  flight. 

A  brief  review  of  the  hypersonic  aerodynamics,  chemical  physics, 
plasma  physics,  and  electromagnetic  wave  scattering  phenomena 
involved  in  the  study  of  plasma  effects  on  radar  cross  sections  is 
given.  It  is  shown  that  the  effect  of  the  plasma  sheath  on  the  radar 
cross  section  of  a  hypersonic  body  can  be  considered  separately 
from  the  effect  of  the  plasma  trail  left  behind  the  body.  Several 
theoretical  models  of  electromagnetic  wave  scattering  from  plasma* 
covered  bodies  are  formulated  and  analyzed.  The  theoretical  results 
of  these  analyses  show  both  absorption  and  enhancement  effects  that 
could  be  caused  by  the  plasma  sheath. 

Theoretical  predictions  of  plasma  sheath  effects  calculated  from  the 
analyses  of  the  theoretical  models  are  compared  with  experimental 
results.  It  is  seen  that  the  measured  absorption  and  enhancement  are 
both  significantly  greater  than  predicted  by  the  theoretical  results. 
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Possible  reasons  for  these  discrepancies  between  theory  and  experi¬ 
ment  are  discussed. 

It  is  concluded  that  the  physical  mechanisms  that  cause  the  absorption 
and  enhancement  effects  observed  experimentally  are  not  yet  com¬ 
pletely  understood  and  consequently  not  completely  included  in  the 
theoretical  models.  Several  possible  additional  factors  are  suggested 
for  incorporation  into  more  advanced  theoretical  models  of  electro¬ 
magnetic  wave  scattering  from  plasma  sheathed  bodies  in  order  to 
predict  the  radar  absorption  effect  more  accurately. 
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SECTION  I 
INTRODUCTION 


One  aspect  of  the  general  study  >f  the  interaction  between  electro¬ 
magnetic  waves  and  plasmas  that  has  received  significant  attention 

recently  is  the  effect  of  a  plasma  sheath  surrounding  a  body  on  the 

(1  through  5)* 

radar  cross  section  of  the  body,  &  Such  a  plasma  sheath 

can  be  produced  around  a  body  in  hypersonic  flight  through  a  gaseous 
atmosphere.  Thus,  interest  in  this  subject  arises  in  the  consideration 
of  the  dynamic  radar  cross  sections  of  reentry  objects.  Since  certain 
characteristic  effects  are  caused  by  the  plasma  sheath,  it  is  also  of 
interest  to  consider  these  effects  as  the  keys  to  new  plasma  diagnostic 
techniques.  For  these  two  reasons,  a  thorough  understanding  of  the 
causative  mechanisms  behind  these  effects  and  an  adequate  theoretical 
analysis  that  can  quantitatively  predict  these  effects  are  necessary. 

In  this  report,  the  particular  phenomenon  known  as  the  radar  absorp¬ 
tion  effect  is  discussed.  The  usual  manifestation  of  this  effect  is  the 
drastic  reduction  of  the  radar  cross  section  of  a  body  when  it  is 
surrounded  by  a  plasma  sheath  having  certain  spatial  configuration 
and  electromagnetic  properties.  This  effect  has  been  observed  during 
the  atmospheric  reentry  of  hypersonic  test  vehicles  v  '  and  in  labora¬ 
tory  studies  of  the  flow  fields  around  bodies  in  hypersonic  flight 
through  various  ambient  gases.  v  h  '  Up  to  the  present  time  a 
quantitatively  adequate  theoretical  analysis  that  predicts  this  effect 
under  all  the  conditions  that  have  been  observed  has  not  been  available, 

although  certain  aspects  of  this  effect  are  v  ^11  understood  and  can  be 

(5) 

predicted  with  reasonably  good  accuracy.  ’ 

♦lumbers  in  parentheses  refer  to  references  listed  at  the  end  of  this 
report. 


1/2 


TR63-217A 


SECTION  II 

DESCRIPTION  OF  PHENOMENA 


Study  of  the  plasma  sheath  effect  on  radar  cross  sections  involves 
interplay  between  the  disciplines  of  hypersonic  aerodynamics,  chem¬ 
ical  physics,  plasma  physics,  and  electromagnetic  wave  scattering. 
The  simultaneous  application  of  the  first  three  disciplines  is  re¬ 
quired  in  order  to  determine  the  configuration  and  physical  properties 
of  the  plasma  sheath  formed  around  a  body  in  hypersonic  flight  through 
an  ambient  atmosphere.  The  physical  properties  of  the  plasma  sheath 
are  in  general  described  by  the  phase  space  distribution  functions  of 
all  the  constituents,  which  simplifies  to  the  specification  of  composi¬ 
tion,  temperature,  and  density  in  the  special  case  of  thermodynamic 
and  chemical  equilibrium.  Important  parameters  entering  into  this 
aspect  of  the  phenomenon  are  the  ambient  gas  density,  temperature, 
and  composition,  and  the  body  configuration,  angle  of  attack,  and 
speed.  The  latter  two  disciplines  are  then  applied  to  determine  the 
electromagnetic  response  of  the  plasma-sheathed  body  to  the  applied 
electromagnetic  wave  from  the  radar  and  thus  determine  the  electro¬ 
magnetic  wave  scattering.  Important  parameters  entering  into  this 
aspect  of  the  phenomenon  are  the  configuration  and  physical  proper¬ 
ties  of  the  body  and  the  surrounding  plasma  sheath,  the  frequency  of 
the  incident  electromagnetic  wave,  and  the  relative  orientation 
between  the  body  and  the  direction  of  propagation  of  the  incident  and 
scattered  electromagnetic  waves.  As  indicated,  the  aerodynamic, 
chemical,  and  plasma  aspects  of  the  problem  can  be  decoupled  from 
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the  plasma  and  electromagnetic  wave  scattering  aspects  in  most  cases, 
in  particular  when  the  energy  density  of  the  electromagnetic  field 
is  negligible  compared  to  the  thermal  and  kinetic  energy  density  of 
the  flow  field.  This  allows  independent  determination  of  the  pro¬ 
perties  of  the  flow  field  and  independent  study  of  electromagnetic 
wave-plasma  interaction  and  scattering  from  plasma-coated  bodies. 
However,  the  final  prediction  of  the  details  of  plasma  sheath  effects 
in  any  given  aerodynamic  situation  necessarily  involves  the  combina¬ 
tion  of  the  results  of  computations  from  all  these  disciplines.  It  is 
important  to  emphasize  that  if  any  one  of  these  computations  is  in 
error,  then  the  final  prediction  will  be  incorrect. 

Consider  a  body  traveling  through  an  ambient  gaseous  atmosphere 
at  hypersonic  speed.  Energy  lost  by  the  body  through  aerodynamic 
forces  acting  on  the  body  is  transferred  to  the  gas  that  is  disturbed 
by  the  body.  This  energy  transfer  is  manifested  in  the  formation  of  a 
strong  shock  wave  around  the  front  of  the  high-speed  body.  The  gas 
in  the  shock  layer  around  the  front  of  the  body  has  received  a  large 
amount  of  energy  that  is  distributed  among  the  various  degrees  of 
freedom  of  the  gas  particles.  Some  excitation,  dissociation,  and 
ionization  of  the  gas  particles  will  result,  thus  producing  radiation 
and  free  electrons.  A  significant  amount  of  the  radiation,  from  the 
infrared  through  the  ultraviolet,  is  radiated  out  of  the  shock  layer. 

The  ultraviolet  radiation  may  produce  additional  ionization  ahead 
of  the  shock  front  (precursor  ionization)  when  it  is  absorbed  by  the 
ambient  gas  particles.  Also,  some  of  the  high-energy  free  electrons 
may  diffuse  ahead  of  the  shock  front,  causing  additional  precursor 
ionization.  As  the  partially  ionized  gas  in  the  shock  layer  flows  from 
the  front  of  the  body  (the  stagnation  region)  around  to  the  sides  of  the 
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body  it  expands  and  cools,  thereby  reducing  the  free  electron  concen¬ 
tration.  This  thin  layer  of  partially  ionized  gas  (plasma)  surrounding 
the  hypersonic  body  constitutes  the  plasma  sheath.  When  the  flowing 
plasma  reaches  the  rear  of  the  body  it  rapidly  expands  to  fill  the 
volume  swept  out  behind  the  body,  and  further  expansion  to  ambient 
pressure  also  occurs  very  rapidly.  The  effect  of  this  expansion  is  to 
reduce  the  free  electron  concentration  in  the  plasma  left  behind  the 
body  (the  trail  or  wake)  very  rapidly*  so  that  shortly  behind  the  body 
this  quantity  may  be  orders  of  magnitude  smaller  than  in  the  stagna¬ 
tion  region.  Farther  behind  the  body  the  free  electron  concentration 
in  the  trail  decreases  by  diffusion,  attachment,  and  recombination 
processes  at  a  much  slower  rate. 

This  simple  picture  of  the  flow  field  may  be  complicated  by  the  form¬ 
ation  of  a  recompression  zone  behind  the  body,  and  the  addition  of 
foreign  particles  to  the  flow  field  by  the  ablation  of  material  from 
the  body  surface.  Determination  of  the  exact  variation  of  the  free 
electron  concentration  as  well  as  other  particle  specie  concentrations 
around  the  body  and  in  the  trail  requires  detailed  consideration,  in 
both  the  equilibrium  and  non-equilibrium  regimes,  of  the  aerodynamic 
and  chemical  processes  involved.  These  considerations  lie  at  current 
frontiers  of  knowledge  in  both  sciences. 

When  an  electromagnetic  wave  propagates  in  a  region  containing  free 
electrons  (such  as  a  plasma)  the  electric  and  magnetic  fields  of  the 
wave  influence  the  motion  of  the  free  electrons,  which  in  turn  extract 
energy  from  the  electromagnetic  wave.  This  absorbed  energy  is 
either  returned  to  the  electromagnetic  field  as  coherent  scattered 
radiation  (with  a  broadened  frequency  spectrum  due  to  doppler  and 
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other  effects)  or  converted  to  thermal  energy  (both  non-coherent 
radiant  energy  and  random  kinetic  energy).  Thermalization  of  the 
energy  is  the  result  of  "collisions"  between  the  free  electrons  and 
other  particles  in  the  plasma.  When  the  free  electron  concentration 
in  the  plasma  is  very  small  the  coherent  scattered  radiation  does 
not  change  the  total  electromagnetic  field  significantly,  the  energy 
losses  are  small  ,  and  an  incident  electromagnetic  wave  propagates 
almost  unaffected  by  the  "transparent"  plasma.  If  the  free  electron 
concentration  is  sufficiently  large,  the  coherent  scattered  radiation 
and  energy  losses  contribute  significantly  to  the  description  of  the 
total  electromagnetic  field. 

In  order  to  facilitate  quantitative  discussion  of  electromagnetic 
wave-plasma  interaction  it  is  convenient  to  introduce  the  concepts 
of  plasma  frequency  and  electron  collision  frequency.  The  plasma 
frequency  is  a  measure  of  the  free  electron  concentration,  and  is 
given  in  Appendix  A-l.  The  physical  significance  of  the  plasma  fre¬ 
quency  is  that,  if  there  are  no  electron  collisons  with  other  particles, 
an  electromagnetic  wave  with  a  frequency  higher  than  the  plasma 
frequency  will  propagate  in  the  plasma,  while  if  the  wave  frequency 
is  lower  than  the  plasma  frequency  propagation  is  not  possible  (the 
electric  and  magnetic  fields  merely  "diffuse"  in  the  plasma).  The 
electron  collision  frequency  is  just  the  average  number  of  collisions 
per  second  between  each  free  electron  and  other  particles  in  the 
plasma  required  to  account  for  the  energy  losses.  From  a  detailed 
viewpoint,  the  concept  of  electron  collision  frequency  is  really 
somewhat  of  an  artificiality.  In  a  plasma  in  which  the  collision 
frequency  is  not  zero  (losses  are  present),  electromagnetic  wave 
propagation  is  always  possible  although  with  greatly  increased 
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attenuation  when  the  wave  frequency  is  less  than  the  plasma  frequency. 
Thus,  if  the  electromagnetic  wave  frequency  is  much  greater  than  both 
the  plasma  frequency  and  the  electron  collision  frequency  the  plasma 
appears  essentially  transparent  to  the  wave,  while  if  the  plasma  fre¬ 
quency  approaches  or  exceeds  the  wave  frequency  the  wave  suffers 
significant  scattering  (and  absorption  if  the  electron  collision  fre¬ 
quency  also  approaches  or  exceeds  the  wave  frequency). 

Now  consider  an  electromagnetic  wave  (such  as  a  radar  pulse)  incident 
on  the  plasma-sheathed  body  and  its  plasma  trail.  If  the  radar  frequency 
is  sufficiently  higher  than  the  stagnation  region  plasma  frequency  and 
electron  collision  frequency,  the  incident  radar  wave  will  be  unaffected 
by  any  of  the  plasma  since  at  all  other  locations  around  and  behind 
the  body  the  plasma  frequency  and  electron  collision  frequency  are 
usually  lower  than  in  the  stagnation  region.  Scattering  of  this  incident 
electromagnetic  wave  will  thus  be  caused  by  the  body"  alone,  and  the 
radar  cross  section  will  be  the  free-space  radar  cross  section  of  the 
body.  If  the  radar  frequency  is  lowered  (or  the  electron  concentration 
is  increased,  due  to  higher  ambient  gas  density  and/or  higher  body 
speed)  the  radar  wave  will  first  be  affected  by  the  plasma  in  the 
stagnation  region,  and  with  continued  radar  frequency  reduction  (or 
continued  electron  concentration  increase)  the  radar  wave  will  be 
affected  progressively  by  more  of  the  plasma  sheath  around  the  sides 
of  the  body.  Throughout  this  regime  the  trail,  although  present,  will 
usually  be  essentially  transparent  to  the  radar  wave  and  consequently 
will  not  contribute  to  the  radar  cross  section;  the  plasma  sheath 
surrounding  the  body  will  then  determine,  to  a  great  extent,  the 
effective  radar  cross  section. 
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If  the  electron  collision  frequency  is  very  small  compared  to  both  the 
radar  frequency  and  plasma  frequency,  the  plasma  sheath  will  appear 
reflective,  approaching  the  appearance  of  a  solid  reflective  coating 
when  the  plasma  frequency  exceeds  the  radar  frequency.  If,  on  the 
other  hand,  the  electron  collision  frequency  is  comparable  to  the 
radar  frequency  and  plasma  frequency,  the  plasma  sheath  will  appear 
absorptive.  Hence,  in  this  regime  the  plasma- sheathed  body  will 
appear  to  be  a  reflective  or  absorptive  object  having  the  configuration 
of  the  outer  surface  of  the  plasma  sheath,  the  expansion  zone  behind 
the  body  may  give  a  small  additional  contribution,  and  the  trail  will 
be  essentially  transparent. 

Further  reduction  of  the  radar  frequency  (or  increase  of  the  electron 
concentration)  will  cause  the  interaction  between  the  incident  electro¬ 
magnetic  wave  and  the  plasma  in  the  trail  to  become  significant,  and 
thus  the  trail  will  begin  to  contribute  to  the  total  radar  cross  section. 
The  effective  length  of  the  trail  in  which  significant  microwave- 
plasma  interaction  occurs  will  tend  to  increase  as  the  radar  frequency 
is  reduced  or  the  electron  concentration  left  behind  the  body  increases. 
It  is  important  to  note,  however,  that  this  interaction  may  take  the  form 
of  either  scattering  or  absorption,  depending  on  the  magnitude  of  the 
electron  collision  frequency,  and  therefore  it  does  not  immediately 
follow  that  the  radar  cross  section  of  the  trail  will  necessarily 
increase  under  these  conditions.  In  the  parameter  regime  in  which 
the  trail  does  contribute  significantly  to  the  total  effective  radar 
cross  section  the  variations  due  to  the  contribution  from  the  plasma- 
sheathed  body  may  be  negligible  in  comparison  to  the  variations  due  to 
the  contribution  from  the  plasma  trail. 
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Consideration  of  the  phenomena  described  in  the  preceding  paragraphs 
suggests  an  initial  subdivision  of  the  problem  into  essentially  indepen¬ 
dent  determinations  of  the  plasma  sheath  effect  on  the  radar  cross 
section  of  the  body  and  the  radar  cross  section  of  the  plasma  trail 
as  an  additional  target  element,  from  which  the  total  effective  radar 
cross  section  of  the  complex  target  can  be  determined  by  appropriate 
summation.  This  subdivision  is  particularly  useful  if  the  significant 
parameter  sensitivities  of  the  sheath  effect  and  the  trail  effect  occur 
in  different  parameter  regimes. 

Determination  of  the  plasma  sheath  effect  on  the  radar  cross  section 
of  the  body  (and  the  radar  cross  section  of  the  plasma  trail)  is  accom¬ 
plished  by  first  the  determination  of  the  physical  properties  of  the 
plasma  and  their  spatial  variation  (the  "boundaries”  of  the  plasma), 
and  then  the  determination  of  the  electromagnetic  wave  scattering 
characteristics  from  which  the  radar  cross  section  is  derived. 

Effective  decoupling  of  the  physical  properties  and  "boundaries" 
determination  from  the  electromagnetic  scattering  determination 
depends  upon  the  condition  that  the  incident  electromagnetic  wave 
does  not  significantly  alter  the  physical  properties  or  the  spatial 
variation  of  the  plasma.  This  is  usually  the  case  in  normal  radar 
applications,  since  the  coherent  energy  supplied  by  the  incident 
radar  wave  is  negligible  compared  to  the  normal  thermal  and  kinetic 
energy  level  in  the  plasma. 

If  the  details  of  the  frequency  spectrum  broadening  of  the  scattered 
coherent  electromagnetic  energy  are  not  to  be  considered,  the 
physical  properties  of  the  plasma  and  the  details  of  the  energy  transfer 
mechanisms  can  usually  be  subsumed  in  a  phenomenological  description 
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of  the  propagation  of  the  incident  electromagnetic  wave.  This 
phenomenological  approach  consists  of  the  assignment  of  effective 
constitutive  parameters  (dielectric  constant  and  conductivity)  to  the 
plasma,  as  well  as  any  other  media  involved,  and  the  subsequent 
treatment  of  the  electromagnetic  behavior  on  a  classical  electro¬ 
magnetic  field  theory  basis.  Thus,  the  plasma  is  considered  to  be  a 
lossy  (and  in  general,  nonlinear,  anisotropic,  and  inhomogeneous) 
dielectric  material. 

The  analytical  determination  of  the  scattering  of  an  incident  electro¬ 
magnetic  wave  by  localized  spatial  distributions  of  material 
characterized  by  electromagnetic  constitutive  parameters  is  a 
classical  boundary-value  problem  in  electromagnetic  field  theory. 

The  complete  description  of  the  scattering  effects  can  be  expressed 
in  terms  of  boundary- matched  sets  of  solutions  to  Maxwell's  equations 
(one  set  of  solutions  for  each  of  the  regions  involved).  Information 
about  the  nature  of  this  solution  in  the  high-frequency  regime  can  be 
obtained  from  a  geometrical  optics  solution  to  the  problem.  The  results 
of  a  geometrical  optics  solution  to  the  problem  must  be  used  with 
caution, however,  since  incomplete  or  highly  misleading  interpretations 
are  possible  if  the  applicability  criteria  are  not  satisfied.  Even  though 
its  applicability  may  be  limited,  a  geometrical  optics  solution  should 
not  be  excluded  a  priori  since  it  can  be  obtained  with  relative  ease 
compared  to  the  rigorous  solution  of  the  electromagnetic  boundary- 
value  problem.  Often,  certain  characteristic  features  of  the  solution 
can  be  easily  seen  from  the  geometrical  optics  solution  even  when  the 
exact  solution  is  extremely  complicated. 
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SECTION  III 

MODELS  OF  PHENOMENA,  ANALYTIC  APPROACHES, 
AND  THEORETICAL  RESULTS 


It  is  obvious  from  the  preceding  discussion  that  the  complete 
physical  situation  is  too  complicated  to  analyze  in  comprehensive 
detail  initially.  Part  of  the  difficulty  of  a  comprehensive  analysis 
lies  in  the  cross-coupling  between  the  different  scientific  disciplines 
involved,  each  of  which  has  its  own  unsolved  problems.  Thus,  to 
provide  an  initial  approach  to  the  determination  of  the  plasma  sheath 
effect  on  the  radar  cross  section  of  a  body,  simple  models 
characterizing  the  main  features  of  the  structure  and  properties  of 
a  plasma-sheathed  body  were  evolved.  Furthermore,  in  order  to 
expedite  the  initial  determination  of  the  behavior  of  these  models  and 
to  simplify  the  interpretation  of  the  analytical  results,  simple  analytic 
approaches  to  the  various  aspects  of  the  behavior  were  used. 

Basic  Model 

The  body  confib  ^ation  chosen  to  represent  typical  blunt-nosed 
hypersonic  vehicles,  for  frontal  aspects,  is  the  sphere.  In  order  to  • 
make  the  electromagnetic  wave  scattering  analysis  tractable  the 
sphere  is  assumed  to  be  perfectly  reflective.  This  represents  a 
realistic  approximation  for  metallic  bodies.  This  assumption, 
however,  makes  the  model  not  directly  applicable  to  the  consideration 
of  ablative -dielectric-coated  vehicles. 
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It  is  assumed  that  all  ionization  around  the  front  of  the  sphere  is 
confined  within  the  shock  layer.  Thus,  the  bow  shock  wave  around 
the  front  of  the  sphere  defines  the  outer  boundary  of  the  plasma  sheath. 
In  this  case,  the  effect  of  precursor  ionization,  if  any,  is  not  con¬ 
sidered. 

The  determination  of  the  exact  shape  of  the  bow  shock  wave  formed 
around  the  front  of  a  body  in  hypersonic  flight  through  a  gaseous 
atmosphere  involves  lengthy  and  complicated  aerodynamic  com¬ 
putation,  even  for  relatively  simple  body  configurations.  In  order  to 
avoid  this  complication,  it  is  assumed  that  the  plasma  sheath  forms 
a  constant-thickness  layer  around  the  sphere.  If  only  the  frontal  view 
of  the  sphere  is  considered  (the  stagnation  region)  this  approximation 
appears  to  be  reasonably  in  accord  with  the  actual  physical  picture. 

(9,  10,  11)  An  explicit,  exact  formula  for  the  calculation  of  the 
plasma  sheath  thickness  at  the  stagnation  point  (shock  wave  detach¬ 
ment  distance)  has  not  been  obtained  up  to  the  present  time;  however, 

(12  13  14) 

approximate  formulas  are  available.  ’  ’ 

The  physical  properties  of  the  plasma  sheath  that  determine  its 

electromagnetic  wave  interaction  behavior  are  subsumed  in  the 

specification  of  an  effective  dielectric  constant  and  conductivity  of 

the  plasma.  Thus,  the  spectral  broadening  of  the  scattered  electro- 
(15) 

magnetic  energy  ’  is  not  considered.  The  effective  electromagnetic 
constitutive  parameters  (dielectric  constant  and  conductivity)  of  the 
plasma  are  derived  from  the  consideration  of  the  propagation  of  a 
small -amplitude  plane  transverse  electromagnetic  wave  through  an 
unbounded  uniform  plasma  with  no  externally  applied  magnetic  field, 
considering  the  effect  of  only  the  free  electrons.  This  analysis  yields 
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linear  and  isotropic  constitutive  parameters,  given  in  Appendix  A-l, 
dependent  only  on  the  plasma  frequency,  electron  collision  frequency, 
and  incident  wave  frequency.  Thus,  the  anisotropic  effects  in  the 
plasma  introduced  by  the  application  of  a  magnetic  field  are  not 
considered.  The  spatial  variation  (homogeneity)  of  the  plasma  con¬ 
stitutive  parameters  throughout  the  plasma  sheath  must  be  known 
explicitly  in  order  to  determine  the  electromagnetic  wave  scattering 
behavior.  This  requires  determination  of  the  spatial  variation  of 
electron  concentration  and  electron  collision  frequency  throughout 
the  plasma  sheath,  which  is  a  complicated  and  lengthy  aerophysical 
computation.  In  order  to  avoid  this  complication,  it  is  assumed  that 
the  constitutive  parameters  of  the  plasma  are  spatially  uniform 
(homogeneous)  throughout  the  plasma  sheath. 

As  a  result  of  the  simplifying  approximations  detailed  in  the  preceding 
paragraphs,  the  initial  theoretical  investigation  of  the  plasma  sheath 
effect  on  the  radar  cross  section  of  a  hypersonic  body  is  reduced  to 
the  problem  of  determining  how  the  radar  cross  section  of  a  perfectly 
reflective  body  is  changed  when  it  is  coated  with  a  uniform-thickness 
layer  of  linear,  isotropic,  homogeneous  lossy  dielectric  material. 

In  the  case  of  the  spherical  body,  an  exact  analytic  solution 
is  available  (see  Appendix  A- 2),  but  the  interpretation  of  the 
mathematical  results  is  very  difficult  without  extensive  numerical 
computation.  In  keeping  with  the  precept  of  a  simplified  approach 
to  the  initial  study  of  the  behavior  of  the  plasma  sheath  effect,  the 
short-wavelength  (or  high-frequency)  regime  was  selected,  thus 
allowing  the  application  of  a  geometrical  optics  approach  to  the 
determination  of  the  radar  cross  section  of  both  the  bare  sphere  and 
the  dielectric -coated  (plasma- sheathed)  sphere.  The  steps  in  this 
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analytic  approach  are  detailed  in  Appendixes  A-3  through  A-7.  This 
approach  considers  the  propagation  of  an  electromagnetic  wave  from 
a  radar  toward  the  dielectric-coated  metallic  sphere  where  partial 
reflection  of  the  wave  back  toward  the  radar  occurs  at  the  outer 
surface  of  the  dielectric  coating.  Partial  transmission  (penetration) 
of  the  incident  wave  into  the  dielectric  coating  also  occurs,  and  this 
transmitted  wave  is  attenuated  as  it  propagates  through  the  dielectric 
coating  to  the  sphere  surface  where  it  is  reflected  so  that  it  propagates, 
with  additional  attenuation,  back  to  the  outer  surface  of  the  dielectric 
coating.  At  the  outer  surface  of  the  dielectric  coating  the  attenuated 
wave  is  partially  transmitted  out  of  the  dielectric  coating  toward  the 
radar  and  partially  reflected  back  into  the  dielectric  coating,  where 
the  attenuated  propagation,  sphere  surface  reflection,  etc.  cycle  is 
repeated  over  and  over  again.  The  total  electromagnetic  wave  re¬ 
turned  to  the  radar  from  the  dielectric-coated  sphere  thus  consists 
of  the  coherent  superposition  of  the  wave  reflection  from  the  dielectric 
surface  and  an  infinite  number  of  decreasing-intensity  wave  contri¬ 
butions  emerging  from  the  dielectric  surface  due  to  the  multiple 
reflections  within  the  dielectric  layer.  The  intensity  of  the  total 
returned  electromagnetic  wave  is  the  measure  of  the  radar  cross 
section.  The  geometrical  optics  approach  used  here  accounts  for 
electromagnetic  wave  interference  and  thus  predicts  scintillation 
effects  that  arise  from  this  phenomenon.  This  geometrical  optics 
approach  is  valid  when  the  size  and  radii  of  curvature  of  the  body 
are  much  larger  than  the  wavelength  of  the  incident  electromagnetic 
wave;  in  general,  it  gives  the  direct  specular  contribution  to  the  radar 
cross  section,  but  does  not  account  for  surface  or  ducted  waves 
which  may  also  contribute  to  the  radar  cross  section. 
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When  this  analysis  is  applied  to  the  situation  where  the  dielectric 
represents  a  plasma  (see  Appendix  A-8),  several  interesting  results 
are  obtained.  In  particular,  the  most  significant  result  of  the 
analysis  is  the  prediction  of  conditions  under  which  the  radar  cross 
section  of  the  plasma-covered  metallic  sphere  vanishes.  These 
conditions  therefore  define  complete  ’’radar  absorption",  and  a  radar 
absorption  effect  will  occur  whenever  these  conditions  are  approached. 
The  amount  of  absorption  (decrease  of  radar  cross  section)  that 
occurs  in  any  given  situation  is  dependent  on  how  closely  the  com¬ 
plete  absorption  conditions  are  approached  in  that  situation.  It  can 
be  seen  from  Appendix  A-8  that  the  normalized  plasma  frequency 
required  for  complete  absorption  is  of  the  order  of  unity  for  small 
(less  than  unity)  values  of  normalized  electron  collision  frequency, 
and  varies  approximately  as  the  square  root  of  the  normalized  electron 
collision  frequency  for  large  (greater  than  unity)  values  of  normalized 
electron  collision  frequency.  The  minimum  plasma-layer  thickness 
for  which  complete  absorption  can  occur  is  of  the  order  of  one-third 
of  the  free-space  wavelength  of  the  radar  wave  when  the  normalized 
electron  collision  frequency  is  large,  arid  increases  at  a  moderate 
rate  as  the  normalized  electron  collision  frequency  is  decreased 
below  unity.  These  results  are  comparable  to  the  results  obtained 
from  the  rigorous  analysis  of  the  case  of  an  infinite, flat,  metallic 
sheet  covered  with  a  uniform  layer  of  lossy  plasma  (see  Appendix 
A-9). 


15 


TR63-217A 


Advanced  Model  - 

In  the  basic  model  considered  in  the  previous  section,  the  spatial 
variation  of  the  plasma  properties  throughout  the  plasma  sheath 
around  the  hypersonic  body  was  neglected.  This  appears  to  be  one  of 
the  most  obvious  deficiencies  of  the  basic  model.  Actually,  the 
plasma  properties  vary  spatially  in  both  the  angular  (azimuthal) 
direction  around  the  body  and  the  radial  direction  away  from  the  body 
surface.  In  general,  both  the  plasma  frequency  and  the  electron 
collision  frequency  decrease  in  the  azimuthal  direction  away  from  the 
stagnation  point.  The  radial  variation  can,  in  general,  be  more 
complicated.  It  is  possible  for  the  plasma  frequency  to  increase, 
reach  a  maximum,  and  then  decrease  in  the  radial  direction  away 
from  the  body  surface.  The  electron  collision  frequency  around 
the  front  of  the  body  is  usually  essentially  constant  from  the  body 
surface  out  to  the  shock  front,  where  it  decreases  abruptly  to  the 
value  for  the  ambient  gas. 

If  chemical  equilibrium  is  maintained  in  the  shock  layer,  the  plasma 
frequency  in  the  shock  layer  around  the  front  of  the  body  also  is 
relatively  constant  from  the  body  surface  out  to  the  shock  front,  where 
it  decreases  rapidly  to  zero  unless  significant  precursor  ionization 
is  present.  Such  precursor  ionization  effectively  extends  the  radial 
range  of  non-zero  plasma  frequency  outside  the  shock  front,  although 
usually  with  a  much  smaller  value  than  in  the  shock  layer  and  with  a 
smoothly  decreasing  magnitude  outward  in  the  radial  direction.  If 
chemical  equilibrium  is  not  maintained  in  the  shock  layer,  two 
different  forms  of  radial  variation  of  the  plasma  frequency  are 
possible.  In  one  case,  the  plasma  frequency  begins  to  increase 
inside  the  shock  front  and  continues  monotonically  increasing 
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toward  the  body  surface.  The  gradient  of  plasma  frequency  depends 
on  the  chemical  reaction  rates  involved  in  the  production  of  free 
electrons  from  the  neutral  gas  particles.  In  the  other  case,  the 
plasma  frequency  also  begins  to  increase  inside  the  shock  front 
but  reaches  a  maximum  from  which  it  decreases  in  the  direction 
toward  the  body  surface.  This  "overshoot"  may  be  significantly 
large  compared  to  the  magnitude  of  the  plasma  frequency  near  the 
body  surface.  The  details  of  this  spatial  variation  again  depend  on 
the  chemical  reactions  occurring  in  the  gas  behind  the  shock  front. 
Here  again,  precursor  ionization  may  be  present  with  a  generally 
similar  effect  to  the  case  where  chemical  equilibrium  is  maintained. 

Some  general  idea  of  how  radial  and  azimuthal  gradients  in  the 
electromagnetic  properties  of  the  plasma  sheath  may  modify  the 
solution  obtained  for  the  basic  model  (uniform  plasma  layer)  can 
be  obtained  by  separate  consideration  of  gradients  in  the  two 
directions.  Here  all  the  simplifications  of  the  basic  model  are 
employed  except  that  the  plasma  layer  is  assumed  to  have  either 

radially  or  azimuthally  non-uniform  electromagnetic  properties. 

(16) 

Some  work  has  been  reported  v  '  on  the  scattering  of  electro¬ 
magnetic  waves  by  a  dielectric-coated  metallic  sphere,  where  the 
dielectric  coating  is  uniform  in  thickness  but  has  radial  variation 
of  its  electromagnetic  properties.  Considering  this  model  in  regard 
to  the  radar  absorption  effect,  it  is  possible  for  a  lossy  dielectric 
layer  with  a  radial  gradient  to  reduce  the  magnitude  of  the  back- 
scattered  electromagnetic  wave,  and  hence  the  radar  cross  section, 
more  than  it  would  be  reduced  by  an  equivalent  uniform  lossy 
dielectric  layer.  In  terms  of  the  plasma  parameters,  if  the  plasma 
frequency  decreases  smoothly  in  the  radial  direction  away  from 
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the  sphere  surface  from  a  value  larger  to  a  value  significantly 
smaller  than  the  radar  frequency  and  the  electron  collision  frequency 
is  moderately  large  compared  to  the  radar  frequency,  considerable 
absorption  will  occur  provided  also  that  the  distance  over  which  the 
plasma  frequency  decreases  is  significantly  large  compared  to  the 
radar  wavelength.  If  the  change  in  plasma  frequency  occurs  over  a 
distance  that  is  a  small  fraction  of  the  radar  wavelength,  no  signi¬ 
ficant  difference  between  this  case  and  the  uniform-plasma-layer 
situation  will  exist.  Thus,  for  a  plasma  sheath  that  is  only  a  small 
fraction  of  a  radar  wavelength  in  thickness,  the  effect  of  a  radial 
gradient  may  be  neglected  by  simply  using  an  average  value  of  the 
plasma  frequency  in  the  layer  for  the  basic  model  analysis. 

The  effect  of  azimuthal  variation  of  the  electromagnetic  properties 
of  the  dielectric  coating  around  a  metallic  sphere  on  the  electro¬ 
magnetic  wave  scattering  characteristics  of  the  dielectric-coated 
sphere  is  analyzed  in  Appendixes  A-10 through  A-13.  In  this  analysis 
all  the  assumptions  made  in  the  analysis  of  the  basic  model  are 
made  except  uniformity  of  the  dielectric  properties,  the  short - 
wavelength  (or  high-frequency)  regime  is  similarly  considered, 
and  the  geometrical  optics  approach  is  also  used  to  determine  the 
radar  cross  section  of  the  non-uniform  lossy  dielectric-coated 
metallic  sphere.  The  variation  of  the  dielectric  properties  is 
assumed  to  be  symmetrical  about  an  axis  of  symmetry, and  the  radar 
cross  section  looking  along  this  axis  is  found.  This  corresponds  to 
looking  at  a  plasma-sheathed  hypersonic  sphere  from  the  nose-on 
aspect. 
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Several  interesting  results  are  obtained  from  this  analysis  of  the 

electromagnetic  wave  scattering  characteristics  of  this  non-uniform 

lossy  dielectric-coated  metallic  sphere.  Here  again,  as  in  the  situation 

where  the  dielectric  layer  is  uniform,  conditions  can  be  found  for  which 

the  radar  cross  section  will  vanish.  These  conditions  are  similar  to 

those  obtained  for  the  uniform  dielectric  situation.  In  addition,  the 

dielectric  layer  acts  as  a  lens  which  can  either  focus  or  de-focus  the 

back-scattered  electromagnetic  wave  and  thus  greatly  increase  as  well 

as  decrease  the  radar  cross  section.  This  is  a  new  effect  not  predicted 

by  the  basic  model  (uniform  dielectric).  The  focusing  mechanism 

potentially  can  explain  increases  in  the  radar  cross  section  which  also 

(7) 

have  been  observed.  Unfortunately,  when  strong  focusing  is  pre¬ 
sent  the  applicability  criteria  of  the  geometrical  optics  approach  are 
violated  and  the  quantitative  results  become  questionable.  This  does 
not  mitigate  the  possibility  for  the  physical  mechanism  (gradients) 
to  produce  such  increases  in  radar  cross  section;  it  simply  means  that 
a  more  refined  analysis  must  be  used  to  obtain  accurate  quantitative 
predictions  in  this  situation. 

In  order  to  determine  appropriate  spatial  variations  of  dielectric 
properties  that  correspond  to  the  variations  of  the  plasma  properties 
in  the  plasma  sheath  around  the  front  of  a  hypersonic  sphere,  some  flow 
field  studies  and  computations  are  necessary.  These  are  outlined  in 
Appendix  A- 14.  The  results  of  these  calculations  show,  as  expected, 
that  the  plasma  frequency  and  the  electron  collision  frequency  in  the 
plasma  sheath  both  decrease  monotonically  in  going  from  the  stag¬ 
nation  point  around  to  the  sides  of  the  body.  However,  the  electro¬ 
magnetic  wave  relative  phase  factor  in  the  plasma  layer  can  either 
decrease  or  increase  initially  in  going  from  the  stagnation  point 
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around  to  the  sides  of  the  body  (see  Appendix  A- 14).  It  is  this 
angular  variation  of  the  relative  phase  factor  that  can  cause  focusing 
or  de-focusing  of  the  back- scattered  electromagnetic  wave  at  nose- 
on  aspect. 

In  general,  when  the  plasma  frequency  is  less  than  the  radar  frequency 
the  angular  variation  of  the  relative  phase  factor  is  such  as  to  cause 
de-focusing,  and  when  the  plasma  frequency  is  greater  than  the  radar 
frequency  the  angular  variation  of  the  relative  phase  factor  is  such  as 
to  cause  focusing.  The  exact  ratio  of  plasma  frequency  to  radar  frequency 
at  which  transition  from  de-focusing  to  focusing  occurs  depends  on  the 
gradient  of  electron  collision  frequency.  Specifically,  when  there  is 
no  angular  variation  of  the  electron  collision  frequency  this  transition 
occurs  when  the  plasma  frequency  is  the  square  root  of  two  times  the 
radar  frequency.  As  the  gradient  of  electron  collision  frequency  in¬ 
creases,  the  ratio  of  plasma  frequency  to  radar  frequency  at  which  this 
transition  occurs  increases.  Since  the  plasma  frequency  in  the  plasma 
sheath  increases  with  increasing  speed  of  the  body,  it  follows  that  the 
de-focusing  effect  will  occur  at  a  lower  speed  than  that  for  which  the 
focusing  effect  will  occur. 

When  the  sphere  flow  field  analysis  and  the  non-uniform  lossy  dielectric- 
coated  metallic  sphere  wave  scattering  analysis  are  combined  (see 
Appendix  A- 15)  and  applied  to  the  plasma- sheathed  sphere  problem, 
several  interesting  results  additional  to  those  obtained  for  the  case 
of  the  uniform  plasma  coating  are  obtained.,  It  is  found  that  the  usual 
flow  field  spatial  variations  of  the  plasma  sheath  electromagnetic  pro¬ 
perties  can  produce  gradients  in  the  equivalent  dielectric  coating  such 
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that  both  de-focusing  and  focusing  of  the  nose-on  back- scattered 
electromagnetic  wave  are  possible „  The  relative  parameter  regimes 
in  which  these  effects  occur  have  been  outlined  in  the  preceeding 
paragraph.  The  net  effect  of  de-focusing  and  focusing  on  the  radar 
cross  section  depends  very  sensitively  on  the  exact  phase  and 
amplitude  relationships  between  Jill  the  component  back- scattered 
waves  that  add  coherently  to  give  the  total  back- scattered  wave 
whose  amplitude  is  the  measure  of  the  radar  cross  section.  Potentially, 
both  de-focusing  and  focusing  gradients  can  either  decrease  or 
increase  the  radar  cross  section  compared  to  that  obtained  for  the 
equivalent  "zero-gradient"  situation.  In  general,  however,  a  de- 
focusing  gradient  cannot  increase  the  radar  cross  section  to  a 
value  greater  than  that  obtained  when  the  plasma  coating  is  perfectly 
reflective.  It  can,  of  course,  greatly  decrease  the  radar  cross 
section  and  consequently  enhance  any  radar  absorption  effect.  A 
focusing  gradient  can  greatly  increase  the  radar  cross  section  above 
the  value  obtained  when  the  plasma  coating  is  perfectly  reflective. 

It  can  also  cause  a  pronounced  decrease  when  the  amplitude  and 
relative  phase  relations  between  the  component  back-scattered  waves 
are  just  right.  Whether  de-focusing  or  focusing  will  have  a  significant 
effect  depends  on  how  much  of  the  net  back- scattered  power  has 
passed  through  the  plasma  layer  where  it  is  affected  by  the  de-focusing 
or  focusing  gradients. 

In  other  words,  if  the  incident  wave  does  not  penetrate  into  the 
plasma  layer  significantly,  it  cannot  be  significantly  affected  by  the 
gradients,  and  consequently  the  back- scattered  wave  cannot  exhibit 
significant  de-focusing  or  focusing  effects.  At  low  body  speeds  the 
plasma  frequency  in  the  plasma  sheath  is  low  and  consequently  much 
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of  the  incident  wave  power  passes  through  the  plasma  layer.  It  is 
under  these  conditions  that  the  usual  radar  absorption  effect  occurs  and 
that  de-focusing  gradients  exist.  Thus,  it  may  be  expected  that  de- 
focusing  gradients  will  usually  enhance  the  absorption  effect.  At  high 
body  speeds  the  plasma  frequency  in  the  plasma  sheath  is  high  and 
consequently  very  little  of  the  incident  wave  power  penetrates  into 
the  plasma  layer.  It  is  under  this  condition  that  focusing  gradients 
exist.  Consequently,  because  there  is  so  little  power  available  to  be 
focused,  it  is  not  likely  that  the  radar  cross  section  will  be  greatly 
increased  by  focusing  gradients. 
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SECTION  IV 

APPLICATION  OF  THEORETICAL  RESULTS  TO 
HYPERSONIC  BODIES  AND  COMPARISON  WITH 
EXPERIMENTAL  RESULTS 

Some  experimental  measurements  of  the  radar  cross  sections  of 
hypersonic  bodies  in  a  free-flight  ballistic  range  have  been  report- 

(7) 

ed,  under  ,  nditions  where  a  radar  absorption  effect  would  be 
expected.  These  observations  confirmed  the  existence  of  the  absorp¬ 
tion  effect  under  controlled  environmental  conditions.  It  was  found, 
however,  that  the  measured  absorption  in  some  cases  was  considerably 
greater  than  predicted  by  the  simplest  theoretical  model,  as  exempli¬ 
fied  by  the  ' 'basic  model"  of  the  previous  section  of  this  work. 
Comparisons  of  the  experimental  results  from  Reference  7  with  the 
theoretical  predictions  from  the  "basic  model"  analysis  of  the 
phenomenon  are  shown  in  Figures  1,  2  and  3. 

The  details  of  the  experimental  program,  whose  results  are  shown 
in  Figures  1,  2  and  3,  are  given  in  full  in  Reference  7  and  will  not 
be  repeated  here.  It  should  be  noted,  however,  that  two  different 
hypersonic  body  configurations  were  used  in  the  experimental  work. 

One  model  (copper -capped  model)  was  a  short  plastic  cylinder  (10  mm 
long  and  20  mm  in  diameter)  whose  forward  end  was  rounded.  This 
rounded  nose  was  covered  with  a  copper  cap  (13  mm  nose  radius) 
which  served  as  a  heat  shield.  The  other  body  (copper  sphere)  was  a 
plastic  sphere  covered  with  a  copper  skin  (15  mm  in  diameter)  which 
also  served  as  a  heat  shield.  The  flight  stability  of  the  copper-capped 
models  was  very  good,  as  evidenced  by  spark  photographs  of  the 
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Figure  1  Nose-On  Radar  Cross  Section  of  a  Copper-Capped  Model 
in  Flight  vs  Velocity,  35  Gc 
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Figure  2  Nose-On  Radar  Cross  Section  of  a  Copper-Capped  Model 
in  Flight  vs  Velocity,  70  Gc 
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models  taken  along  the  flight  path.  The  structural  integrity  of  the 
copper  heat  shields  of  both  types  of  models  throughout  their  flights 
was  verified  at  the  end  of  each  flight  by  flash  X-ray  photographs.  In 
the  execution  of  the  experimental  program  the  ambient  air  pressure 
in  the  flight  chamber  was  held  constant  at  10  mm  Hg,  The  models 
were  observed  in  flight  by  a  CW  doppler  radar  at  the  velocities 
shown  by  the  data  points.  The  vertical  extent  covered  by  each  data 
point  shows  the  extent  of  the  variation  of  the  measured  radar  cross 
section  over  the  length  of  the  flight  path  under  observation 
(approximately  30  feet). 

An  outline  of  the  method  of  calculation  of  the  theoretical  predictions 
shown  in  Figures  1,  2  and  3  is  as  follows.  In  order  to  use  the 
"basic  model"  analysis  to  calculate  the  nose-on  radar  cross  section 
of  a  plasma- sheathed  blunt  metallic  body,  it  is  necessary  to  know  the 
plasma  frequency  and  the  electron  collision  frequency  at  the  stagna¬ 
tion  point  and  the  thickness  of  the  plasma  layer  around  the  nose  of 

(17) 

the  body.  Charts  are  available  '  showing  the  flow  field  properties 
of  a  blunt  body  in  hypersonic  flight  through  air  over  a  range  of 
ambient  air  density  and  flight  velocity,  and  these  charts  were  used  to 
obtain  the  plasma  frequency  and  electron  collision  frequency  at  the 
stagnation  point  as  functions  of  body  velocity  for  an  ambient  air 
pressure  of  10  mm  Hg.  The  derivation  of  the  charts  is  based  on  the 
condition  that  chemical  equilibrium  is  maintained  in  the  flow  field. 

In  order  to  check  the  applicability  of  these  charts  for  the  small 
bodies  used  in  the  experimental  work,  special  non-equilibrium  flow 
field  calculations  were  carried  out  for  these  small  body  sizes  at 
specific  velocities  within  the  range  encountered  in  the  experimental 
work.  It  was  found  (see  Appendix  A- 14)  that  near-equilibrium 
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conditions  prevailed  in  the  flow  field  near  the  body  around  the  nose 
region  of  the  bodies,  and  consequently  the  chart  values  of  plasma 
properties  are  reasonably  accurate  even  for  the  small  bodies  used 
in  the  experimental  work.  The  plasma  frequency  and  electron  collision 
frequency  at  the  stagnation  point  as  functions  of  body  velocity  at  an 
ambient  air  pressure  of  10  mm  Hg,  obtained  from  the  charts  and 
used  in  the  theoretical  calculations,  are  shown  in  Appendix  A-14. 

The  value  of  plasma  layer  thickness  used  in  the  calculations  is  the 
actual  measured  value  of  shock  detachment  distance  as  a  function  of 
body  velocity,  taken  from  Reference  7  and  shown  in  Appendix  A-14; 
this  value  is  in  good  agreement  with  other  published  theoretical  and 
experimental  work 5  ^ 

Using  the  above  input  data,  the  theoretical  prediction  curves  of 
nose-on  radar  cross  section  as  a  function  of  body  velocity  shown  in 
Figures  1,  2  and  3  were  calculated  from  equation  (8. 1)  of  Appendix  A-8, 
which  is  the  geometrical  optics  approximation.  This  calculation 
embodies  the  "basic  model"  assumption  that  the  plasma  layer  can  be 
considered  to  be  spatially  uniform.  It  can  be  seen  that  in  all  cases 
the  measured  absorption  is  greater  than  predicted  by  this  theoretical 
approach. 

In  order  to  remove  any  doubt  about  the  validity  of  the  use  of  the  geo¬ 
metrical  optics  approximation  in  deriving  the  theoretical  curves 
shown  in  Figures  1,  2  and  3,  the  rigorous  field  theory  solution  of  the 
"basic  model"  was  obtained  (see  Appendix  A-2).  Using  the  same  input 
data  as  for  the  geometrical  optics  calculation,  the  variation  of  the 
nose -on  radar  cross  section  as  a  function  of  body  velocity  was  calcu¬ 
lated  from  the  rigorous  field  theory  solution,  and  is  shown  in  Figures 
4  and  5. 
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Figure  4  Wave  Solution  for  Copper-Capped  Model,  35  and  70  Gc 
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Figure  5  Wave  Solution  for  Copper  Sphere,  35  and  70  Gc 
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It  can  be  seen  from  Figures  4  and  5  that  there  is  little  significant 
difference  between  the  rigorous  field  theory  approach  and  the  geo¬ 
metrical  optics  approach  in  theoretically  predicting  the  amount  of 
absorption  according  to  the  "basic  model"  assumptions.  It  appears, 
therefore,  that  the  "basic  model"  is  deficient  in  characterizing 
the  essential  mechanisms  that  cause  the  large  amount  of  absorption 
in  the  parameter  regime  covered  by  the  experimental  results. 

It  was  pointed  out  earlier  that  one  of  the  obvious  deficiencies  in  the 
"basic  model"  is  that  it  does  not  account  for  the  gradients  in  the 
plasma  properties  in  the  plasma  sheath  around  the  nose  of  the  body. 
Such  gradients  could,  in  principle,  cause  focusing  or  de-focusing 
of  the  back- scattered  radar  wave.  The  "advanced  model"  was  for¬ 
mulated  to  incorporate  this  additional  mechanism,  and  the  nose-on 
radar  cross  section  of  this  model  was  found  using  the  geometrical 
optics  approach  (see  Appendix  A-15).  In  order  to  use  this  analysis 
to  predict  the  amount  of  absorption  that  will  occur,  it  is  necessary 
to  have  analytic  expressions  for  the  gradients  of  plasma  frequency 
and  electron  collision  frequency  in  the  flow  field  in  the  angular 
direction  away  from  the  stagnation  point  in  the  vicinity  of  the  stag¬ 
nation  point.  Some  discussion  and  analysis  of  this  aspect  of  the 
situation  is  given  in  Appendix  A- 14.  It  is  shown  there  that  the 
variations  of  plasma  frequency  and  electron  collision  frequency  in 
the  angular  direction  near  the  stagnation  point  can  be  approximated 
by  parabolic  functions  of  the  angle  around  the  nose  away  from  the 
stagnation  point.  Using  the  coefficients  of  these  parabolic  represen¬ 
tations  obtained  from  the  non-equilibrium  flow  field  analysis  (see 
Appendix  A- 14)  together  with  the  plasma  properties  at  the  stagnation 
point  used  in  the  previous  calculations,  the  variation  of  the  nose-on 
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radar  cross  section  as  a  function  of  body  velocity  predicted  by  the 
’’advanced  model"  was  calculated  from  Equation  (15. 1)  of  Appendix 
A- 15  and  is  shown  in  Figures  6  and  7. 

It  can  be  seen,  by  comparison  between  Figures  1,  2  and  3  and 
Figures  6  and  7,  that  the  "advanced  model"  predicts  a  somewhat 
different  absorption  behavior  as  a  function  of  sphere  velocity,  but  it 
does  not  account  for  the  very  large  amount  of  absorption  observed 
experimentally.  Some  enhancement  effects  due  to  focusing  of  the  back- 
scattered  radar  wave  are  evident  in  the  "advanced  model"  theoretical 
results.  It  should  be  noted  that  these  theoretical  results  are  based 
on  the  geometrical  optics  approach.  In  view  of  the  good  agreement 
between  the  results  of  the  geometrical  optics  approach  and  the 
results  of  the  rigorous  wave  theory  approach  in  the  analysis  of  the 
"basic  model",  it  is  reasonable  to  expect  the  geometrical  optics 
approach  to  be  applicable  in  the  analysis  of  the  "advanced  model" 
also.  Since  the  "advanced  model"  results  do  not  predict  a  large 
amount  of  absorption  it  appears,  therefore,  that  the  presence  of 
angular  gradients  in  the  plasma  properties  is  not  the  significant 
mechanism  causing  the  large  amount  of  absorption  in  the  parameter 
regime  covered  by  the  experimental  results. 


32 


Figure  6  Geometrical  Optics  Solution  for  Copper-Capped  Model 
35  &  70  Gc 
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Figure  7  Geometrical  Optics  Solution  for  Copper  Sphere,  35  & 
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SECTION  V 

DISCUSSION  OF  RESULTS 


The  comparisons  between  the  experimental  results  available  from 
ballistic  range  studies  of  the  effect  of  the  plasma  sheath  on  the  nose- 
on  radar  cross  sections  of  blunt  metallic  bodies  and  the  theoretical 
predictions  based  on  the  theoretical  models  of  the  phenomenon  for¬ 
mulated  in  this  report  show  that  distinct  discrepancies  between  theory 
and  experiment  exist.  Anomalous  absorption  and  enhancement  are 
both  seen  in  the  experimental  results.  That  is,  the  maximum  absorption 
measured  is  significantly  greater  than  predicted  by  either  the  "basic 
model"  or  the  "advanced  model",  and  the  maximum  enhancement  in 
the  cases  where  it  appears  is  significantly  greater  than  predicted 
by  the  "advanced  model".  These  results  seem  to  indicate  that  either 
the  numerical  values  of  the  plasma  parameters  used  in  the  theoretical 
calculations  were  in  error  or  that  the  theoretical  models  are  deficient 
in  representing  a  dominant  mechanism  involved  in  the  phenomenon. 

Such  deficiencies  may  involve  either  (or  both)  the  aerodynamic  or 
(and)  the  electromagnetic  aspects  of  the  situation. 

Consider  first  the  effects  of  variation  of  the  numerical  values  of  the 
plasma  parameters  used  in  the  theoretical  calculations.  The  para¬ 
meters  involved  are  the  plasma  frequency,  the  electron  collision 
frequency,  and  the  thickness  of  the  plasma  layer.  The  effects  of 
altered  plasma  properties  on  the  theoretical  predictions  of  the 
"basic  model"  are  shown  in  Figures  8,  9  and  10.  The  plasma 
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Figure  9  Effect  of  Plasma  Frequency  and  Collision  Frequency 

Changes  on  Radar  Absorption,  Copper-Capped  Model,  70  Gc 
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frequency  and  electron  collision  frequency  have  been  increased  in 
such  proportions  as  to  increase  the  maximum  absorption  at  approxi¬ 
mately  a  fixed  velocity.  Factors  of  increase  in  plasma  frequency  and 
electron  collision  frequency  of  up  to  10  and  100,  respectively,  are 
considered.  It  is  apparent  from  Figures  8,  9  and  10  that  this  alone 
is  not  sufficient  to  resolve  the  experimental-theoretical  discrepancy, 
and  such  large  deviations  from  the  values  used  in  the  original  calcu¬ 
lations  are  not  at  all  justified.  It  appears,  therefore,  that  the  dis¬ 
crepancy  is  not  due  to  the  use  of  incorrect  values  of  plasma  frequency 
or  electron  collision  frequency  in  the  theoretical  calculations.  The 
effects  of  variation  of  the  numerical  value  of  the  plasma  layer  thick¬ 
ness  on  the  theoretical  predictions  of  the  ,rbasic  model”  are  shown  in 
Figures  11,  12  and  13.  The  value  of  thickness  used  in  the  calculations 
shown  is  increased  by  factors  of  2,  3  and  4  above  the  values  used  in 
the  original  calculations.  These  curves  show  that  maximum  absorption 
of  the  same  order  of  magnitude  as  the  experimental  results  is  pre¬ 
dicted  when  the  plasma  sheath  thickness  used  in  the  calculations  is 
approximately  three  times  the  value  originally  used.  This  is  consistent 
with  the  requirements  for  vanishing  radar  cross  section  given  in 
Appendix  A-8,  since  the  tripled  thickness  is  approximately  one  or 
two  halves  of  the  radar  wavelength  at  35  Gc  or  70  Gc,  respectively, 
for  the  copper-capped  model  and  approximately  one-half  the  radar 
wavelength  at  70  Gc  for  the  copper  sphere.  It  should  be  noted  that 
the  numerical  values  of  plasma  sheath  thickness  used  in  the  original 
theoretical  prediction  calculations  are  actually  the  shock  detachment 
distances  measured  from  photographs  of  the  actual  flow  fields.  These 
measured  shock  detachment  distances  are  in  good  agreement  with 
other  published  data  and  cannot  be  in  error  by  more  than  20  percent. 
Hence,  if  the  ionization  is  confined  to  the  region  behind  the  shock 
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Figure  11  Effect  of  Changes  in  Stand-Off  Distance  on  Radar  Absorption, 
Copper-Capped  Model,  35  Gc 
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Figure  12  Effect  of  Changes  in  Stand-Off  Distance  on  Radar  Absorption, 
Copper-Capped  Model,  70  Gc 
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Figure  13  Effect  of  Changes  in  Stand-Off  Distance  on  Radar 
Absorption,  Copper  Sphere,  70  Gc 
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front,  the  plasma  layer  is  too  thin  (by  a  factor  of  approximately 
three)  to  account  for  the  measured  absorption. 

The  theoretical  models  formulated  here  of  the  phenomena  involved 
in  the  effect  of  the  plasma  sheath  on  the  radar  cross  sections  of 
blunt  metallic  bodies  have  not  considered  the  possibility  of  the 
presence  of  ionization  outside  of  the  shock  front  around  the  nose  of 
the  body.  Such  precursor  ionization  could  be  produced  by  either 
(or  both)  electron  diffusion  ahead  of  the  shock  front  or  (and)  ultra¬ 
violet  photoionization  of  the  air  ahead  of  the  luminous  shock  front. 
Precursor  ionization  ahead  of  the  shock  front  effectively  increases 
the  plasma  sheath  thickness  to  a  value  greater  than  the  shock  de¬ 
tachment  distance.  The  plasma  frequency  in  a  precursor  ionization 
region  caused  by  either  electron  diffusion  or  by  photoionization  will 
decrease  smoothly  in  the  radial  direction  away  from  the  body.  Such 
a  transition  region  between  free  space  and  the  shock  front  around 
the  nose  of  the  body  can,  if  the  parameters  are  appropriate, 
significantly  influence  the  scattering  of  an  incident  electromagnetic 
wave  and  consequently  affect  the  radar  cross  section. 

Two  mechanisms  by  which  such  a  transition  region  can  significantly 
decrease  the  radar  cross  section  are  (1)  the  actual  absorption  of  the 
incident  wave  power  in  the  transition  region  (real  absorption)  and 
(2)  the  reflection  and  refraction  or  diffraction  of  the  incident  wave 
in  the  transition  region  into  directions  other  than  the  back-scattering 
direction  (deviative  absorption).  Real  absorption  requires  that  the 
transition  region  be  lossy,  which  implies  that  the  electron  collision 
frequency  as  well  as  the  plasma  frequency  are  not  negligible  compared 
to  the  radar  frequency.  Some  analysis  of  this  situation  was  given  in 
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Reference  7>  where  it  was  shown  that  in  order  to  produce  significant 
absorption  in  a  transition  region  with  constant  electron  collision  fre¬ 
quency  and  exponentially  decreasing  plasma  frequency  the  e-folding 
distance  had  to  be  greater  than  one -fourth  of  the  radar  wavelength, 
and  the  peak  plasma  frequency  and  electron  collision  frequency  had  to 
be  greater  than  one -tenth  of  the  radar  frequency. 

An  analysis  of  the  deviative  absorption  situation  has  not  yet  been  made, 
but  it  is  reasonable  to  expect  in  this  situation  also  that  in  order  to 
produce  a  significant  effect  the  peak  plasma  frequency  wouldhaveto 
be  at  least  one-tenth  of  the  radar  frequency  or  greater.  A  more  refined 
examination  of  this  matter  is  planned.  At  the  present  time  there  is  a 
serious  lack  of  extensive  quantitative  information  about  the  amount  of 
precursor  ionization  that  may  be  present  at  various  flight  conditions. 
However,  some  preliminary  measurements  of  precursor  ionization  have 
been  made  v  '  and  it  appears  from  these  measurements  that  the  peak 
plasma  frequency  in  the  prec  ursor  region  is  too  low  at  the  flight 
conditions  encountered  in  the  experimental  work  to  give  significant 
interaction  effects.  At  higher  velocities  and  for  lower  radar  fre¬ 
quencies  the  effect  of  precursor  ionization  will  be  more  significant, 

and  may  account  for  some  of  the  absorption  effects  observed  during 

(6) 

the  atmospheric  reentry  of  hypersonic  test  vehicles  '  . 

Because  the  plasma  layer  is  so  thin  compared  to  the  wavelength  of 
the  incident  radar  wave  in  the  experimental  regime  studied  in 
Reference  7  ,  gradients  within  the  layer  in  the  radial  direction  away 
from  the  body  could  not  greatly  influence  the  overall  scattering 
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behavior.  Hence,  the  fact  that  such  gradients  have  not  been  included 
in  the  theoretical  models  could  not  greatly  alter  the  conclusions 
drawn  so  far.  In  principle,  however,  such  gradients  could  considerably 
affect  the  scattering  behavior  when  they  extend  over  a  region  which 
is  of  the  order  of  a  wavelength  in  size. 

Another  physical  mechanism,  not  incorporated  in  any  of  the  theoretical 
models  so  far,  that  could  possibly  account  for  both  the  anomalous 
absorption  and  the  increases  in  radar  cross  section  is  the  coupling 
of  the  transverse  electromagnetic  radar  wave  with  longitudinal 
plasma  waves  (electroacoustic  waves).  The  anomalous  absorption 
could  result  when  strong  coupling  between  the  waves  exists  and  the 
electroacoustic  waves  are  highly  damped,  thus  forming  an  effective 
’’sink"  for  the  coherent  energy.  Increases  in  the  radar  cross  section 
could  result  if  the  electroacoustic  waves  resonate  within  the  plasma 
sheath.  These  ideas  are  speculative  at  this  time,  but  work  is  planned 
to  substantiate  (or  deny)  these  possibilities. 

The  possibility  is  under  study  that  peculiar  diffraction  effects  may 
cause  the  anomalous  absorption  and  the  increase  in  radar  cross 
section.  It  is  possible  that  the  copper-capped  plastic  model  may 
give  rise  to  some  strong  real  diffraction  effects  both  with  and  without 
the  plasma  sheath  present.  Studies  of  this  model  and  configurational 
variations  of  it  are  being  made  in  a  static  radar  cross  section  range 
in  order  to  answer  this  question. 
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SECTION  VI 
CONCLUSION 


Several  theoretical  models  of  radar  wave  scattering  from  plasma- 
sheathed  blunt  metallic  bodies  have  been  formulated  and  analyzed. 

The  results  of  the  analyses  show  that  the  theoretical  models  potentially 
can  account  for  both  extreme  absorption  and  enhancement  under 
appropriate  conditions.  Predictions  of  the  effects  of  the  plasma  sheath 
on  radar  cross  sections  based  on  these  theoretical  models  have  been 
compared  with  experimental  results.  The  theoretical  models  fail  to 
predict  correctly  the  amount  of  absorption  and  enhancement  observed 
experimentally  for  the  particular  values  of  parameters  involved.  The 
apparent  reason  for  failure  of  the  theory  to  give  numerical  results 
consistent  with  experimental  results  is  that  the  plasma  layer  thickness 
used  in  the  numerical  calculations,  as  deduced  from  the  most  reasonable 
aerodynamic  considerations,  is  too  small  to  yield  significant  absorption 
or  enhancement  effects.  It  appears,  therefore,  that  the  physical 
mechanism  that  causes  the  absorption  and  enhancement  effects  observed 
experimentally  is  not  yet  completely  understood  and  consequently  is 
not  included  in  the  theoretical  models.  Further  work  on  the  theoretical 
models  is  necessary  to  incorporate  into  them  the  additional  factors 
mentioned  earlier  in  this  work,  and  experiments  over  a  wider  range 
of  parameters  should  be  carried  out  to  better  define  the  theoretical- 
experimental  relationship. 
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.  .  APPENDIX  A-l 

ELECTROMAGNETIC  CONSTITUTIVE 
PARAMETERS  OF  PLASMA 


In  the  phenomenological  viewpoint  of  the  interaction  between 
electromagnetic  waves  and  matter  it  is  customary  to  consider  all 
nonionized  media  to  be  completely  characterized  (electromagneti- 
cally)  by  the  constitutive  parameters 

1.  permeability  (  ), 

2.  capacitivity  (  €  ),  and 

3.  conductivity  (  ). 

It  is  highly  advantageous  to  extend  this  concept  to  include  ionized 
media,  particularly  plasma.  Plasma  is  essentially  a  gaseous  mixture 
of  electrically  neutral  and  charged  particles  of  atomic  or  sub- atomic 
size  exhibiting  essentially  no  net  charge  density.  If  a  plasma  has  a 
sufficiently  large  number  of  each  type  of  particle  per  unit  applicable 
volume,  the  concept  of  configuration  space  density  has  meaning  and 
the  velocity  distribution  of  the  particles  in  a  unit  applicable  volume 
can  be  described  by  a  velocity  space  density.  The  configuration  and 
velocity  space  densities  (phase  space  density)  describe  the  dynamic 
properties  of  the  plasma.  If  the  time  variation  of  the  phase  space 
density  from  causes  other  than  the  incident  electromagnetic  wave  is 
very  small  during  one  cycle  of  the  incident  wave,  the  concept  of 
constitutive  parameters  can  be  applied  to  the  plasma. 
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The  electromagnetic  constitutive  parameters  of  the  plasma  are 
assumed  here  to  be  the  same  as  the  effective  constitutive  parameters 
derived  under  the  condition  of  a  small- signal  plane  transverse  electro 
magnetic  wave  propagating  through  an  unbounded  uniform  (homo¬ 
geneous)  plasma  with  no  externally  applied  electric  or  magnetic 
fields  (isotropic),  considering  the  effects  of  only  the  free  electrons 
and  assuming  a  constant  electron  collision  frequency.  This  derivation 
is  well  documented  in  the  literature/^  the  results  of  which  give  the 
plasma  permeability  /V 

a.  i)  ^ 

where 


jXj  =  free- space  permeability, 


the  plasma  capacitivity 


(1.2) 


€  6. 

\f  cr 


1 


i 

i 


where 


£Kr  =  free-space  capacitivity, 

6r  =  plasma  dielectric  constant, 

Ne  =  plasma  electron  concentration, 
=  electron  charge, 
nnNe  =  electron  mass, 
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<o  =  incident  electromagnetic  wave  angular  frequency 
(2  IT  f), 

v'c  =  plasma  electron  collision  frequency, 
and  the  plasma  conductivity  o'*,  , 


(1.3) 


Ne  ‘te8'  *4 

(  60  *  ♦  t41) 


These  equations  can  be  written  in  more  compact  form  by  introducing 
the  relative  conductivity  , 


(1.4) 


€vr  ^ 


and  the  plasma  frequency  a) 


9 


(1.5) 


Ne 


‘/e 


The  plasma  dielectric  constant  €.  y,  and  relative  conductivity 
(S’  r  can  now  be  expressed  in  the  forms 


(1.6)  €r 


(1.7)  or 

r 


“3  W*  + 
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It  is  convenient  to  introduce  a  normalized  plasma  frequency 
given  by 

(1.8)  rtjtt 

^  CD 


and  a  normalized  electron  collision  frequency  Cl  c  given  by 

(i.9)  n,  * 

to 

in  terms  of  which  equations  (1.  6)  and  (1. 7)  can  be  written  in  the 
forms 


(1.10) 


jlA 

I  -V  il 


z 

C 


(1.11) 


JLfcl. 

i  +  n.1 


(2) 

It  can  be  shown  ’  that  the  electromagnetic  wave  propagation  factor 
for  a  plane  transverse  electromagnetic  wave  propagating 
through  an  unbounded,  homogeneous,  isotropic  plasma  characterized 
by  the  above  three  electromagnetic  constitutive  parameters  is 
given  by 

(1.12)  tp  =  0C+  +  j  (er>^)4(er"iaV),i 
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where  the  attenuation  factor  CL  jz,  is  given  by 


(1.13)  CCm  = 


1 Arl 

vr 


% 


-  €. 


% 


and  the  phase  factor  /Q^  is  given  by 


(1.14)  ^ 


u)  A  v  j  Kl 

VT“ 


% 


The  free-space  wavelength  of  an  electromagnetic  wave  is 


(1.15) 


2  TT  _  2.  IT _ 

A-  "  “!6w>,rl  ^ 


which  can  be  used  to  simplify  equations  (1. 12),  (1. 13),  and  (1. 14)  to 


(1.16)  /  = 
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(1.17)  OC^  = 


vr-fT 

Air 


-er  + 

. 


(1.18)  /$f>  ~ 


ir-sr 

Air 


«r 


It  is  also  convenient  to  define  a  relative  attenuation  factor  OC.^ 
given  by 


(1.19)  oC,. 


VT 


“€r  + 


e,*  + 


S', 


and  a  relative  phase  factor  given  by 


In  terms  of  the  normalized  plasma  frequency  XI  and  the  norma¬ 
lized  electron  collision  frequency  -Qc  ,  equations  (1. 16),  (1.17), 
and  (1. 18)  are 


(i.2i)  =  j  ii 

Ay 


SI Jf 


C 
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so  that  the  relative  attenuation  factor  cCY  and  the  relative  phase 
factor  /8p  for  a  plasma  can  be  written  in  the  forms 


The  relationship  between  the  electromagnetic  wave  wavelength  in 
the  plasma  and  the  free -space  wavelength  Xy.  is 
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(1.26) 


(1)  E.  C.  Jordan,  Electromagnetic  Waves  and  Radiating  Systems, 
Prentice -Hall,  1950,  pp.  659—662. 

(2)  Bendix  Systems  Division,  Wave  Propagation,  Reflection, 
and  Transmission  in  a  Lossy  Semi-Infinite  Cold  Plasma7 
By  G.if.  flesher,  Research  Note  8,  Ann  Arbor,  Michigan, 

20  June  1960. 
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APPENDIX  A-2 

ELECTROMAGNET!  C  WAVE  BACK- SCATTERING  CROSS 
SECTION  OF  A  METALLI C  SPHERE  COATED  WITH 
A  UNIFORM  LOSSY  DIELECTRIC 


The  rigorous  wave  theory  determination  of  the  scattering  of  an 
incident  linearly  polarized  plane  transverse  electromagnetic  wave 
by  a  uniform  isotropic  spherical  body  with  a  concentric  uniform 
isotropic  spherical  coating,  both  of  which  may  have  arbitrary  con¬ 
stitutive  parameters,  immersed  in  an  arbitrary  uniform  isotropic 
surrounding  medium,  has  been  reported  in  the  literature.  ^  This 
solution  can  be  specialized  to  the  case  where  the  inner  spherical 
body  is  a  perfect  conductor  (metal),  the  coating  has  arbitrary 
dielectric  constant  and  conductivity  (lossy  dielectric)  but  a  perme¬ 
ability  the  same  as  free  space,  and  the  surrounding  medium  is  free 
space.  The  back- scattered  plane  transverse  electromagnetic  wave 
is  linearly  polarized  in  the  same  direction  as  the  incident  wave.  The 
electromagnetic  wave  back-scattering  cross  section  (radar  cross 
section)  q*  of  this  special  composite  configuration  is  given  by 


(2.1) 


C f  - 


4  TT 


/_  f 

'ns  i 


-I 


2.  ♦ 


•  I  I 


!  Z. 


^-'Yl  w/rt 


where 


free -space  wavelength  of  the  incident  wave,  and 
scattering  amplitude  coefficients. 


The  scattering  amplitude  coefficients  are  given  by 
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(2.2)  a^  = 


(2.  3)  4 ^ 


o 

Z»OU) 

z43)(nU) 

0 

Z^YNc) 

zf(Mc) 

zY3)  (m<0 

-  n  *U) 

*u  (*u) 

nf(wi) 

o 

'uYNc) 

*1  i  VnO 

CYfU) 

zfU) 

z“(Ma) 

Z^YnJ) 

0 

1-h*  (Nc) 

li^Nc) 

-  N1?U) 

i4°Cna> 

i) 

-  z,<»U) 

z«(iU) 

z^YlU) 

The  factor  N  is  the  complex  index  of  refraction  of  the  lossy  dielectric 
coating,  and  is  given  by 
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where 


*r 

O', 


_£A  _ 

~  €vr 

=  «±- 

G\r  <0 


= 


dielectric  constant  of  dielectric 

relative  conductivity  of  dielectric 

capacitivity  of  dielectric 

conductivity  of  dielectric 

capacitivity  of  free  space 

angular  frequency  of  incident  wave  (  £  It  f  ). 


The  quantities  c  and  d  in  the  arguments  of  the  functions  in  equation 
(2. 1)  are  defined  by 

(2.5)  c  =  a  ^  r 

A  \r 


(2.6) 


2  IT  (  r  +  ft.) 

A.  r 


where 

r  =  radius  of  metal  sphere 

ft.  =  thickness  of  dielectric  coating 


The  function  X.^  (x)  is  the  spherical  Bessel  function  of  the  first 
kind,  of  integer  order  n ,  and  complex  argument  X.  This  function 
is  related  to  the  ordinary  Bessel  function  of  the  first  kind,  of  order 
X  ♦  Vz ,  and  complex  argument  X  by 


6)  v 

(2.7)  X  = 


tr 


\Vk 


sx/ 
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The  function  is  the  spherical  Bessel  function  of  the  third 

kind,  of  integer  order  y\  ,  and  complex  argument  X  .  This  function 

is  related  to  the  ordinary  Hankel  function  of  the  second  kind,  of 

(2) 

order  -nW*  ,  and  complex  argument  X  by  ' 


(2.8) 

2^1  X)  -  ( 

-rr  ^ 

*x\ 

1  tw(x| 

The  functions 

are  defined  by 

(2.9) 

II 

X 

0  £ 
s- 

Ix^'1  lx)]' 

X 

(2.  10) 

,>!  = 

[xi'Jl  lx)]' 

X 

where  the  prime  indicates  differentiation  with  respect  to  the 
argument  X. 

The  complicated  nature  of  this  solution  precludes  observation  of 
the  behavior  of  the  radar  cross  section  as  a  function  of  the  input 
variables  r  ,  k,  €T  ,  ,co,  and  merely  by  inspection. 
Computation  on  an  electronic  digital  computer  is  the  only  practical 
way  in  which  to  obtain  numerical  values  of  the  solution  for  study, 
because  of  the  lengthy  calculation  involved.  Even  with  this  approach, 
there  are  parameter  regimes  in  which  the  solution,  as  given  in  the 
form  of  equation  (2. 1)  et  seq,  does  not  converge  well  and  consequently 
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the  accuracy  of  the  numerical  result  is  poor.  This  occurs  when  r 
is  very  large  compared  to  X  ^  and/or  the  magnitude  of  N  is  very 
small  compared  to  unity.  In  these  parameter  regimes  other  analytic 
forms  of  the  solution  must  be  used. 

In  order  to  adapt  the  above  analytic  solution  to  the  special  case 
where  the  dielectric  coating  represents  a  lossy  plasma,  it  is  merely 
necessary  to  express  N  in  terms  of  the  plasma  variables,  giving 


(2.11) 


.  I  flJ-  '  ** 

Kl  =  I - = -  i 

'  t  '  j  -fic  ' 


where  and  nc  are  defined  in  Appendix  1. 


(1)  A.  L.  Aden  and  M.  Kerker,  "Scattering  of  Electromagnetic 
Waves  from  Two  Concentric  Spheres,"  Journal  of  Applied 
Physics.  Vol.  22,  No.  10,  October  1951,  pp.  1242-1246. 

(2)  A.  Erdelyi,  editor,  Higher  Transcendental  Functions,  Vol,  n, 
McGraw-Hill,  1953,  pp.  9  -  10. 
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APPENDIX  A-3 

GEOMETRICAL  OPTICS  APPROACH  TO  THE 
DETERMINATION  OF  ELECTROMAGNETIC  WAVE 
SCATTERING  CROSS  SECTION 

The  electromagnetic  wave  scattering  cross  section  <f  of  an  object  is 
defined  in  terms  of  the  ratio  of  the  scattered  wave  average  power  den¬ 
sity  at  the  receiver  location  to  the  incident  wave  average  power 
density  at  the  object  location  ”5.  by 

2  Sr 

(3.1)  4  r  4 It  ““ nr 

Si 

where 

V  distance  between  object  location  and  receiver  location. 

In  order  for  this  definition  to  be  complete,  the  polarization  characteris¬ 
tics  of  the  incident  and  scattered  waves  at  the  object  and  receiver  loca¬ 
tions,  respectively,  must  also  be  specified.  When  the  incident  wave  is 
linearly  polarized  in  some  reference  direction,  then  in  general  the 
scattered  wave  will  have  a  component  polarized  in  this  same  direction 
(parallel -polarized)  and  a  component  polarized  in  the  orthogonal 
direction  (cross-polarized). 

The  definition  given  in  equation  (3. 1)  can  be  written  to  give  the  parallel- 
polarized  and  the  cross-polarized  components  of  the  scattering  cross 
section  by  using  the  parallel -polarized  and  the  cross-polarized  compo¬ 
nents  of  the  scattered  wave  average  power  density  at  the  receiver  location. 
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The  analytical  determination  of  the  scattering  of  an  incident  electro¬ 
magnetic  wave  by  a  localized  macroscopic  region  with  known  constitu¬ 
tive  parameters  and  spatial  configuration  is  a  boundary-value  problem 
in  electromagnetic  field  theory.  The  solution  of  the  scattering  problem 
gives  the  scattered  wave  at  a  receiver  location  when  the  incident  wave 
at  the  object  location  is  specified.  The  complete  description  of  the 
scattering  effects  can  be  expressed  in  terms  of  boundary- matched  sets 
of  solutions  to  Maxwell's  equations  (one  set  of  solutions  for  each  of  the 
regions  involved).  Information  about  the  nature  of  this  solution  can  be 
obtained  from  a  geometrical  optics  solution  to  the  problem.  The  results 
of  a  geometrical  optics  solution  to  the  problem  must  be  used  with 
caution, however,  since  incomplete  or  highly  misleading  results  are 
possible  if  the  applicability  criteria  are  not  met.  Even  though  its 
applicability  may  be  limited,  a  geometrical  optics  solution  should  not 
be  excluded  a  priori  since  it  can  be  obtained  with  relative  ease  com¬ 
pared  to  the  rigorous  solution  of  the  electromagnetic  boundary-value 
problem. 

It  can  be  shown  ^  from  Maxwell's  equations  that  the  geometrical 
optics  (ray  theory)  solution  is  reasonably  accurate  when 

(1)  the  fractional  change  in  the  constitutive  parameters  over  a 
distance  of  one  wavelength  is  small  (for  continuously  varying 
constitutive  parameters)  or  if  the  total  change  in  the  constitu¬ 
tive  parameters  takes  place  over  a  distance  very  much  smaller 
than  one  wavelength  (for  an  abrupt  boundary  between  two 
regions) ,  and 

(2)  the  fractional  change  in  the  spacing  between  adjacent  rays 
over  a  distance  along  the  rays  of  one  wavelength  is  small.  This 
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means  that  the  geometrical  optics  solution  is  open  to 
question  if 

(a)  the  rays  pass  through  a  focus,  or 

(b)  the  rays  are  reflected  or  refracted  at  a  surface  with  a 
size  or  radius  of  curvature  that  is  not  large  compared 
to  the  wavelength. 

The  geometrical  optics  approach  can  be  used  to  determine  the  total 
electric  field  intensity  at  a  receiver  location,  from  which  the  average 
power  density  can  be  determined.  To  accomplish  this,  it  is  necessary 
to  consider 

(1)  the  length  of  each  ray  path  from  the  incident  wave  front  to  the 
receiver  location  (to  obtain  the  transmission  phase  shift  and 
attenuation,  and  the  change  in  polarization  characteristics  if 
anisotropic  media  are  traversed), 

(2)  the  magnitude  changes,  phase  shifts,  and  polarization 
characteristics  changes  that  occur  at  the  boundaries  encountered 
by  each  ray  path  (expressed  by  the  reflection  and  transmission 
factors),  and 

(3)  the  magnitude  changes  that  occur  along  each  ray  path  due  to 
geometrical  divergence  of  the  ray  bundle  (expressed  by  the 
geometrical  divergence  factor) 

in  order  to  determine  the  magnitude,  phase,  and  polarization  relation¬ 
ships  between  the  components  that  add  up  to  give  the  total  electric  field 
intensity  at  the  receiver  location. 

Using  the  ray  path  field  intensity  addition  technique,  the  total  scattered 
wave  average  power  density  arriving  at  a  receiver  location  ”5^  is  given 
by 


A3-3 


TR63-217A 


where 

E„  =  peak  value  of  total  electric  field  intensity  at  receiver 
location  (linearly  polarized) 

ZOR  =  characteristic  impedance  of  medium  (lossless)  in  which 
receiver  is  located. 


The  total  electric  field  intensity  at  the  receiver  location  ER  is  the  sum 
of  the  contributions  from  all  the  ray  paths  that  terminate  at  the  re¬ 
ceiver  location. 


(3.3)  Ew  a  £  Er4i 

*fc*» 

where 

ERk  =  electric  field  intensity  of  kth  ray  path  at  receiver  location 
(linearly  polarized) 

Note  that  equations  (3. 2)  and  (3.3)  can  be  written  for  each  of  two 
orthogonally  polarized  electric  field  intensities. 


In  the  scattering  situation,  all  ray  paths  that  reach  the  receiver  location 
originate  from  incident  wave  ray  paths,  hence  the  initial  ray  path 
electric  field  intensity  of  all  ray  paths  is  the  incident  wave  (linearly 
polarized)  electric  field  intensity  E..  The  incident  wave  average  power 
density  is  given  by 


(3.4) 


x 


«Zoi. 
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where 

2  =  characteristic  impedance  of  medium  (lossless)  in  which 

wave  source  is  located. 

In  general,  each  ray  path  will  encounter  boundaries  at  which  it  is 
either  reflected  or  transmitted,  and  consequently  the  magnitude,  phase, 
and  polarization  characteristics  of  the  ray  path  electric  field  intensity 
will  be  affected  by  the  reflection  or  transmission  factors  at  these  bound¬ 
aries.  The  curvature  of  the  boundaries  encountered  by  each  ray  path 
introduces  an  overall  geometrical  divergence  factor,  which  affects  the 
ray  path  average  power  density.  Transmission  along  each  ray  path 
introduces  a  phase  shift  between  the  incident  wave  electric  field  inten¬ 
sity  and  the  electric  field  intensity  of  that  ray  path  at  the  receiver  lo¬ 
cation.  This  phase  shift  can  be  measured  between  two  conveniently 
located  reference  planes,  one  perpendicular  to  the  incident  wave  pro¬ 
pagation  direction  and  in  the  same  medium  as  the  incident  wave  source, 
and  the  other  perpendicular  to  the  scattered  wave  propagation  direction 
and  in  the  same  medium  as  the  receiver  location  but  located  a  fixed 
distance  from  the  receiver  location  (preferably  as  close  to  the  scattering 
body  as  possible).  Transmission  through  lossy  media  causes  attenuation 
of  the  electric  field  intensity  along  the  ray  path.  Transmission  through 
anisotropic  media  causes  changes  in  the  polarization  characteristics 
along  the  ray  path. 

When  all  the  media  involved  are  isotropic,  uniform  (changes  in  con¬ 
stitutive  parameters  occur  only  at  sharp  boundaries  between  regions), 
and  no  changes  in  polarization  characteristics  occur  at  the  boundaries, 
a  general  formulation  for  the  scattering  cross  section  can  be  made. 
Simultaneous  consideration  of  the  effects  of  all  reflections,  transmissions, 
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geometrical  divergence,  phase  shifts,  and  losses  along  the  kin  ray 
path  (in  which  there  are  a  total  of  reflections  at  the  surfaces  r, 
a  total  of  transmissions  at  the  surfaces  t,  a  total  of  uniform 
media  traversed,  a  total  of  A^  uniform  lossy  media  traversed,  and 
the  source  and  receiver  locations  are  in  the  same  lossless  medium), 


gives 


(3.5) 


Ei 


Vz 


* 

rr  <°i 


fcr 


I>iR* 


t'\ 


•j 


T-fe  /w  B* 
TT  ~4t  TT 


t  =  I 


IT 

OL*  I 


-oc 


where 


OC 


fi- 


iRk 

7 

kr 

’kt 

ka 

Jka 

kb 

Jkb 


geometrical  divergence  factor  for  kth  ray  path 
r  surface  reflection  factor  in  the  k  ray  path 
t**1  surface  transmission  factor  in  the  k^  ray  path 
attenuation  factor  along  k**1  ray  path  in  a**1  medium 
length  of  k  ray  path  in  a  medium 

t  Vi  iU 

phase  factor  along  kin  ray  path  in  b  medium 
length  of  kin  ray  path  in  b  medium 


Equations  (3.  2)  through  (3.5),  for  and  equal,  give 


(3.6)  = 

Si. 
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(3.7) 


2 

*  -  Vz 
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TT 

"  “dt 
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TT 

ZL 

4*i 
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Thus  the  electromagnetic  wave  scattering  cross  section  can  be  ex¬ 
pressed  in  terms  of  a  geometrical  optics  solution  of  the  scattering 
problem,  using  a  ray  path  field  intensity  addition  technique. 


(1)  D.  E.  Kerr,  Propagation  of  Short  Radio  Waves.  McGraw-Hill, 
1951,  pp.  53-58. 
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APPENDIX  A-4 

ELECTROMAGNETIC  WAVE  REFLECTION  AND 
TRANSMISSION  ATLOSSV  DIELECTRIC  BOUNDARIES 

It  is  assumed  that  the  boundaries  between  the  lossy  dielectric  and  the 
adjacent  media  are  plane,  and  furthermore  that  all  the  media  involved 
can  be  completely  characterized  (electromagnetically)  by  the  three 
constitutive  parameters:  capacitivity  €  ,  permeability  ,  and  con¬ 
ductivity  c X  .  The  media  bounding  the  lossy  dielectric  are  assumed 
to  be  either  free  space  (a  lossless  dielectric)  or  infinitely  conductive 
material  (perfect  reflector),  and  the  lossy  dielectric  (a  dielectric  with 
nonzero  conductivity)  is  assumed  to  have  free  space  permeability. 

The  reflection  and  transmission  of  electromagnetic  waves  at  plane 
interfaces  between  dielectrics  (both  lossless  and  lossy)  and  perfect 
reflectors  have  been  given  some  attention  in  the  literature.  ^  The 
generalized  Snell's  Law  gives  the  directions  of  the  reflected  and 
transmitted  waves  at  the  interface  when  the  direction  of  the  incident 
wave  is  specified.  The  propagation  direction  vectors  of  the  incident 
wave,  reflected  wave,  transmitted  wave,  and  the  surface  normal  are 
co-planar,  this  common  plane  being  defined  as  the  plane  of  incidence. 
This  is  shown  in  Figure  A4-1  for  a  wave  propagating  from  free  space 
through  a  lossy  dielectric  slab  and  into  free  space  again.  The  gen¬ 
eralized  Snell's  Law  for  the  relationships  shown  in  Figure  A4-1  gives 


TR63-217A 


Figure  A4-1  Reflection  -  Refraction  Geometry  for  Interfaces 
Between  Free  Space  and  Lossy  Dielectric 
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(4.2) 

(4.3)  Aim  9tl  = 

(4'4>  A'  fl 

Aim  0tsf 

where 


Sti  *  flu  *  Of  A 

_ V&~  -Aim.  0{.vr 


6,  +  Aim*  0  +  |  <r/  +  («,-  Aim1  ^  1  ' 

L  ^  j 


S  i  ^  Aim  A.  j 

c  + 

f,  A  +  ' 

Vz.  ' 

V r  u 

r  1 

£  =  dielectric  constant 

tr 

<y*r  =  relative  conductivity. 

The  dielectric  constant  and  relative  conductivity  are  given  by 


(4.5) 

(4.6) 


co 


where 


€  ^  =  dielectric  capacitivity 
€^.  =  free -space  capacitivity 

GO  =  incident  electromagnetic  wave  angular  frequency  (  2ltf  ) 
<Tj  =  dielectric  conductivity. 


The  electromagnetic  wave  reflection  ( /O )  and  transmission  ( ) 
factors  for  an  interface  between  two  media  relate  the  magnitude  and 
phase  of  the  reflected  (Er)  and  transmitted  (Et)  electromagnetic  wave 
electric  field  intensities  at  the  interface  to  the  incident  electromag¬ 
netic  wave  electric  field  intensity  (E^)  at  the  interface,  that  is, 
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(4.7) 

(4.8) 


at  interface 


at  interface 


where  Er,  E^,  and  Ei  are  complex  quantities.  The  reflection  and 
transmission  factors  for  an  interface  are,  in  general,  functions  of 
the  constitutive  parameters  of  the  media  on  the  two  sides  of  the 
interface,  the  angle  of  incidence,  and  the  electric  field  polarization 
with  respect  to  the  plane  of  incidence.  For  normal  incidence 
( 0)  only  the  first  of  the  above  three  functional  dependencies 
enters  the  equations  explicitly.  The  reflection  factor  (  /° &  )  and  the 
transmission  factor  (#t”|fc)  for  an  electromagnetic  wave  propagating 
normally  from  medium  (1)  to  medium  (2)  are  given  by 


(4.9) 

/O  = 

zoz.  ~ 

*oi 

•  1  z 

Zot  * 

Zo. 

(4. 10) 

T..  - 

2  Z  oa. 

u  it  - 

«1 

0 

N 

Z  oi 

where 

Z0,  =  characteristic  impedance  of  medium  (1), 
Zat  =  characteristic  impedance  of  medium  (2). 

The  characteristic  impedance  of  a  lossless  dielectric  is 
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the  characteristic  impedance  of  a  lossy  dielectric  is 


(4. 12) 


- 


*  -i 


<r 

~  i 


*4 


and  the  characteristic  impedance  of  a  perfect  reflector  is 


(4.13)  Ze  -  O 


The  reflection  factor  )  and  transmission  factor  for  an 

electromagnetic  wave  propagating  normally  from  free  space  into  a 
lossy  dielectric  having  free- space  permeability,  using  equations 
(4.9)  through  (4. 12),  are 


(4. 15) 


The  reflection  factor  and  transmission  factor  CTTAv)  f°r  311 

electromagnetic  wave  propagating  normally  from  a  lossy  dielectric 
having  free-space  permeability  into  free  space,  using  equations 
(4. 9)  through  (4. 12),  are 


(4- 16) 
(4. 17) 
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It  can  be  seen  from  equations  (4. 14)  through  (4. 17)  that 


(4. 18) 

/C 

=  ~  f° j-A. 

(4.19) 

Tvri 

=■  1  +  f°  vrdL 

(4. 20) 

"  l  “  t°\ \r<L 

The  reflection  factor  </**.)  and  transmission  factor  <tif)  for  an 
electromagnetic  wave  propagating  normally  from  a  lossy  dielectric 
to  a  perfect  reflector,  using  equations  (4.9),  (4. 10),  (4. 12)  and 
(4.13),  are 


(4.21)  fa  -  _  , 

(4. 22)  li,  *  O 


In  the  special  case  where  the  lossy  dielectric  represents  a  plasma, 
it  is  merely  necessary  to  express  €.r  and  a \  in  terms  of  the 

plasma  variables.  Thus,  the  reflection  factor  (  P^  )  for  an  electro¬ 
magnetic  wave  propagating  normally  from  free  space  into  a  plasma  is 
given  by 


<4-  23>  rt  -  /% 


l  - 

r,  _Qj__ 1 
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l  1  ~  0  j 

\  -*■ 

[  I  -  4  SL  c  j 

•/«, 

A4-6 


TR63-217A 


where  xi*  and  XI  C  are  the  plasma  variables  and  are  defined 
in  Appendix  1.  Equations  (4. 18)  through  (4.  20)  allow  all  the  other 
pertinent  reflection  and  transmission  factors  for  the  interface  to  be 
expressed  in  terms  of  ,  just  as  in  the  case  of  the  lossy 

dielectric. 


(1)  J.  A.  Stratton,  Electromagnetic  Theory.  McGraw-Hill,  1941, 
pp.  500—524. 
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APPENDIX  A-5 

RAY  PATH  GEOMETRY  FOR  BACK-SCATTERING  BY  A 
UNIFORM  LOSSY  DIELECTRIC-COATED  METALLIC  SPHERE 

It  is  assumed  that  the  scattering  object  consists  of  a  perfectly 
reflective  sphere  coated  with  a  uniform  layer  of  lossy  dielectric 
material.  The  ray  path  geometry  of  a  plane  electromagnetic  wave 
incident  on  these  concentric  spherical  boundaries  and  the  resultant 
scattered  wave  ray  path  geometry,  necessary  for  the  determination 
of  the  scattering  cross  section,  can  be  completely  defined  by  the  ray 
paths  in  the  plane  of  incidence  (as  defined  in  Appendix  4).  A  normal 
view  of  the  plane  of  incidence,  showing  the  equatorial  section  of  the 
spherical  boundaries,  is  shown  in  Figure  A5-1.  The  angles  of  refraction 
of  the  ray  paths  shown  are  indicative  of  a  dielectric  constant  with  a 
value  less  than  unity  in  the  coating.  Three  characteristic  ray  path 
configurations  are  apparent,  but  only  the  contributions  to  the  scattered 
wave  due  to  "single-transit”  paths  through  the  dielectric  are  detailed 
(with  the  exception  of  the  "vds^"  path,  which  makes  an  infinite  number 
of  transits  through  the  dielectric).  It  should  be  noted  that  contributions 
to  the  scattered  wave  also  come  from  rays  that  have  been  carried 
around  the  back  of  the  sphere  in  the  dielectric  coating.  The  geometrical 
details  of  each  of  the  three  characteristic  ray  paths  must  be  known 
in  order  to  determine  the  relationship  between  the  incident  wave  and 
the  scattered  wave. 

Consider  the  "v"  ray  path,  in  which  the  ray  is  reflected  from  the 
outer  surface  of  the  dielectric  coating.  This  path  is  detailed  in  Figure 
A5-2. 
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Figure  A5-1  Ray  Path  Geometry  in  the  Plane  of  Incidence 
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-  -  -  Figure  A5-2  Characteristic  Ray  Path  "v" 
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Appendix  4  gives 

(5-1)  =  ©i  \r 

and  the  law  of  sines  for  a  plane  triangle,  with  equation  (5. 1),  gives 

(5.2)  — .  r  .+  *• - .  =  _ 5a _ 

j  2.  $LXr  -  9^  Aim.  |  TT  -  9iv-} 

which  can  be  expanded  and  arranged  to  give 

(5.3)  u *.  ( *.  elw.  -  0ir )  =  W  Qiv. 

The  "far-field”  condition  gives 

(5.4)  Aim  (fc0i,r-0i.r)  =  oo  A**  ©lirj  -  0 

which  means  that 

(5.5)  eif  =  z  9w 

The  length  of  the  "v"  ray  path  is  most  conveniently  measured  from  an 
incident  wave  reference  plane  perpendicular  to  the  direction  of  pro¬ 
pagation  of  the  incident  wave  and  tangent  to  the  outer  surface  of  the 
dielectric  coating,  as  shown  in  Figure  A5-2.  The  lengths  of  the  ray 
path  in  free  space  and  in  the  dielectric  from  this  reference  plane  to 
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the  receiver  location  must  both  be  determined.  From  Figure  A5-2 
it  is  obvious  that  the  length  of  the  ray  path  in  the  dielectric  (LDy) 
is  zero. 


(5.6)  L  Dv  =  0 

The  ray  path  length  in  free  space  from  the  incident  reference  plane 
to  the  dielectric  surface  (  )  is  given  by 

(5.7)  Jfc,  =  (r  +  4}  |  I-  cc*  0lNr) 

and  the  ray  path  length  in  free  space  from  the  dielectric  surface  to 
the  receiver  location  (  Jtt  ),  using  equation  (5. 1),  is  given  by 

(5.8)  C4.  Tr-0-w)  ♦  (£'&*•&,») 

It  is  convenient  to  establish  a  scattered  wave  reference  plane  perpen¬ 
dicular  to  the  direction  of  propagation  of  the  scattered  wave  and  located 
a  distance  from  the  receiver  location.  The  ray  path 

length  from  the  dielectric  surface  to  this  reference  plane  ( )  is 
given  by 

(5.9)  s  (f +  4  j  j  l  *  C**  @i,v)  t  |  Co*  |  2 ®ir!  ~  l 

The  "far-field"  condition  can  be  evaluated  through  the  use  of  equation 
(5.3)  and  1 'Hospital's  rule,  giving 


A5-5 


TR63-217A 


(5.10)  r^»  i*  =  (*  +  •#■)  (l-  0U) 

The  total  ray  path  length  in  free-space  between  the  two  reference 
planes  (LVv)»  for  the  Mfar-field"  condition,  is  given  by 

(5.11)  Lv,r  =  2-  ©L\r) 

For  the  back- scattering  condition  (  Q ^  ~  0  )  it  follows  that 

(5. 12)  LvV  =  0 

Consider  the  "vds"  ray  path,  in  which  the  ray  is  transmitted  through 
the  outer  dielectric  surface,  passes  through  the  dielectric  coating  until 
it  meets  the  sphere  surface,  is  reflected  from  the  sphere  surface, 
again  passes  through  the  dielectric  coating  until  it  meets  the  outer 
dielectric  surface,  and  is  transmitted  through  the  outer  dielectric 
surface  to  the  receiver  location.  In  general,  such  a  ray  path  may 
traverse  the  dielectric  coating  2m  times  before  emerging.  This  path 
is  detailed  in  Figure  A5-3.  Appendix  4  gives 

(5.13)  Aim  0t ^  » 

(5. 14)  = 


_ Msl  -Aa/tv  u  j\r _ 

|  er  +  Airtv*-  eiv  *  p;2  ♦  |e,  -  Qivf  J,/l  |/l- 
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I 


RECEIVER  LOCATION 


Figure  A5-3  Characteristic  Ray  Path  "vds" 
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(5.15)  Uak 


+  1 

c  *  *\  M 

r  ' 

r  ^0u  1  . 

By  symmetry,  it  can  be  seen  that 


(5.16)  «9ti  =  9U 
and 

(5.17)  at,  =  eL, 

The  law  of  sines  gives 


(5.18) 


*•  +  -ft. 


Aim  ♦  0*vr  +  Z  -  ^Lr)  AAm.  (IT-  0^) 


which  can  be  simplified  to 


(5, 19) 


v  + 


RR 


{  Z  <9lvr  +  2  0*  -  0<>)  -  9L[r 


The  "far-field"  condition  gives 


(5.20)  Aim  U  ®t*  +  2  -  »ir)  =  r^,  [  ^  A**  ft,]  =  0 


which  means  that 


(5.21) 


&Lr  s  ^  +  ) 


tv 
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in  the  specific  case  shown  in  Figure  A5-3,  but  in  general 

(5.22)  *  2  (0L<r  +  **.£*.)  nm.  »  1*2*3, 

for  multiple  transits  through  the  dielectric  coating.  The  law  of  sines 
also  gives 


(5.23) 


v  A. 


Aim,  Q. 


tA 


Mj*.  If  -  9t  A  “ 


which  can  be  reduced  to 


(5.24) 


lr+  JL 


Aim  0 


tdL 


A Am.  (  A  + 


and  solved  to  give 

Co4.fi,.  =  [J^l  W1  fiti 


(5.25) 


*  1 1- Pta'14 


i-l^lVVal* 


Using  equation  (5. 13)  in  equation  (5. 25)  gives 


(5.26)  ^  s 


<C,*  +  ( «r  "  ^ 


•  »(■ 

I 


€.y-  Aim4, 0tg.  4* 


fiir + j  <r^  +  fit)rj*  -*  i ■£*&• 


+1  v**1* 
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which  can  be  inverted  to  give  B*.  for  substitution  into  equation 
(5.  22)  to  give  0  iRm  in  terms  of  0^  ,  in  the  form 


(5.27) 


0LR rm  -  z  &i 


\r 


2  rm  ate  c«4.  +  0t(f +|  +  | 

2  f  jLtijw2  eilr  +•[[€„-  9i<r  +[a;1+(€lr-W  9LSY1 
[[e,. + 1»!  6m  *  [tf,.1  +  ( W  9lv ) z  J  Z  j^|2’iwizi9ilr| /z 


The  length  of  the  "vds"  ray  path  is  measured  between  the  same  two 
reference  planes  established  for  the  "v"  ray  path.  From  Figure  A5-3 
it  can  be  seen  that  the  ray  path  length  in  the  dielectric  (  L  0\r<is  )> 
for  one  round  trip,  is  given  by 


(5.28)  ^Dvrds 


Z  r 


Alov  0a. 
Aim.  ©taL 


In  general  the  "vds”  ray  path  makes  m  round  trips  through  the 
dielectric  coating,  thus  the  total  ray  path  length  in  the  dielectric 

(LD»-As<iJisSiven  by 


(5.29) 


L 


D  vds  rm. 


Aim. 

AW  &td. 


which  can  be  expanded,  thru  the  use  of  equations  (5. 13)  and  (5. 26),  to 
give 
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(5.30)  L  OvAsrm  =  2 


|e>  +**'**  9Lv  +  [<*  *  +  (£,.-  lift1  ®i«rf]!4|" 


For  the  back-scattering  condition  (  0  ),  equation  (5.30) 

reduces  to 


<5*31>  L  Dvrd.Srm.  =  nnrL-^- 

The  ray  path  length  in  free  space  from  the  incident  reference  plane  to 
the  dielectric  surface  (  J2|  )  is  given  by 

(5.32)  A,  =  |  r  +  -ft. j  j  l  —  coi.  0iv) 

and  the  ray  path  length  in  free  space  from  the  dielectric  surface  to  the 
receiver  location  (£2)*  using  equation  (5. 17),  is  given  by 

(5.33)  =  (*+•«.)  ccAjir-fli.,,)*  Rrc»4.  (  z9„*  Z&.-0ir] 

The  ray  path  length  from  the  dielectric  surface  to  the  scattered  wave 
reference  plane  (|g),  using  equation  (5.33),  is  given  by 
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(5.34) 


»  r  +  Aj  (  I- 0Lvr )  4 


The  "far-field"  condition  can  be  evaluated  through  the  use  of  equation 
(5.19)  and  1 'Hospital's  rule,  giving 

(5.35)  r^oo  1  (*  +  4-)  (!-«•*■  ®Ur) 

The  total  ray  path  length  in  free  space  between  the  two  reference 
planes  (LVvdgm),  for  the  "far-field”  condition,  is  given  by 

(5.36)  Lwdsm'2  l  r  l|_  ^ 

For  the  back- scattering  condition  (  (9^^^  =  O  ),  equation  (5. 36) 
reduces  to 


<5  37>  LVvdsm  =  0 

Consider  the  "vd"  ray  path,  in  which  the  ray  is  transmitted  through  the 
outer  dielectric  surface,  passes  through  the  dielectric  coating  until  it 
again  meets  the  outer  dielectric  surface,  and  is  transmitted  through  the 
outer  dielectric  surface  to  the  receiver  location.  In  order  for  such  a  ray 
path  to  contribute  to  the  back- scattering  it  must  encircle  the  metal 
sphere  at  least  once.  When  the  dielectric-coating  thickness  is  small 
compared  to  the  sphere  radius,  then  the  ray  path  must  meet  the  outer 
dielectric  surface  many  times  in  going  around  the  sphere.  At  each  of 
these  internal  reflections  the  intensity  is  decreased,  so  that  after  going 
completely  around  the  back  of  the  sphere  the  intensity  is  greatly  reduced 
and  hence  the  contribution  from  such  a  ray  path  may  be  neglected. 
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APPENDIX  A-6 

DIVERGENCE  FACTOR  FOR  BACK- SCATTERING  BY  A 
UNIFORM  LOSSY  DIELECTRIC-COATED  METALLIC  SPHERE 


The  geometrical  divergence  factor  (between  two  points  on  a  ray  path) 
is  defined  as  the  limit  of  the  ratio  of  the  cross  sections  (at  the  two 
points)  of  a  ray  bundle  containing  the  specified  ray  as  the  bundle  shrinks 
to  the  single  specified  ray.  The  physical  significance  of  the  geometrical 
divergence  factor  is  that  it  accounts  for  the  variation  of  the  average 
power  density  along  a  ray  bundle  carrying  constant  power  due  to  the 
variation  of  the  cross- sectional  area  of  the  bundle  along  its  length  for 
the  limiting  condition  of  zero  cross-sectional  area  of  the  bundle  (a 
single  ray)^  If  the  first  point  lies  in  an  initial  plane  perpendicular 
to  the  ray  path  and  the  second  point  lies  in  a  receiver  plane  perpendicular 
to  the  ray  path, the  geometrical  divergence  factor  between  these  two 
locations  (D^)  is  given  by 


(6.1) 


DiR 


4m  A  r, 
Ai  0  A  i 


where  is  the  cross-section  area  of  a  ray  bundle  in  the  initial  plane 
and  A^  is  the  cross-section  area  of  the  same  ray  bundle  in  the  receiver 
plane.  This  divergence  factor  is  independent  of  the  shape  of  the  ray 
bundle  cross  section,  hence  any  form  convenient  to  the  geometry  of  the 
problem  may  be  used. 
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It  is  assumed  here  that  the  scattering  object  consists  of  a  perfectly 
reflective  sphere  coated  with  a  uniform  layer  of  lossy  dielectric 
material.  The  geometrical  divergence  factor  for  a  plane  electro¬ 
magnetic  wave  incident  on  these  concentric  spherical  boundaries, 
necessary  for  the  determination  of  the  scattering  cross  section,  is 
determined  by  the  ray  path  geometry  (as  derived  in  Appendix  5).  The 
geometrical  configuration  used  to  define  a  typical  ray  bundle  that 
can  be  used  in  the  explicit  formulation  of  equation  (6. 1)  for  concentric 
spherical  boundaries  is  shown  in  Figure  A6-1. 

The  initial  area  A.  is  plane  and  parallel  to  the  xy  plane,  hence  it  is 
perpendicular  to  the  incident  rays  from  the  source.  Figure  A6-2 
details  A^,  which  is  the  area  of  the  non- mutual  part  of  a  sector  of  two 
concentric  circles.  It  may  be  seen  from  Figure  A6-1  that 

(0.2)  /O  s  (v+A.)  Aim. 

and 

(6.3)  /°+  A/°  =  (v  +  A)  Aim.  (  9iv  +  A  ) 

from  which  it  follows  that 

(6*4)  Ai.  -  ^  (»•+«&)  +  AS^r)  "Aim  ^iu*J 

The  receiving  area  A^  is  on  the  surface  of  a  sphere  with  center  at  the 
object  location  and  radius  Rp,  hence  it  is  perpendicular  to  the  scattered 
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PLANE  OF  INCIDENCE=  x  z 
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Figure  A6-1 
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rays  from  the  object  when  RR  is  very  much  larger  than  the  dimensions 
of  the  object.  Figure  A6-3  details  A^,  which  is  a  fraction  of  the  area  of 
the  spherical  zone  defined  by  the  planes  z  =  a  and  z  =  b.  It  may  be 
seen  from  Figure  A6-3  that 

(6.5)  a  = 

and 

(6.6)  b=  Ger4.(#-L^ 

from  which  it  follows  that 

(6.7)  Ar  =  A  ^  Rr*  [  ***»  ^LR  “  (  #1*  ♦  A  0|,r  | 

Using  equations  (6.4)  and  (6.7),  the  geometrical  divergence  factor 
for  the  concentric  spherical  boundaries  becomes 

(6*8)  D-  -  Z  Kr1  ^  Cgfc  0LR  -  <**.  (  &iR  +AgvR  1 

|r-»-4)x  A^-d  W  j^+Atf^.)  -x*feLr 

In  which  it  is  necessary  to  express  Si  r,  explicitly  in  terms  of 
Q  tv  before  the  limit  can  be  taken.  For  the  back- scattering  condition, 
equation  (6. 8)  reduces  to 
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Figure  A6-3 


Details  of  Receiving  Area  A 
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(6,9) 


d«.r  = 


y  +  4,)2- 


Xi**.  i  -  ca&  A^LR 
A 9tf~0  AiM,*  A &L\r 


The  two  different  functional  relationships  between  and 

are  given  in  Appendix  5  for  the  two  different  characteristic  ray  paths 
that  contribute  to  the  back- scattering  cross  section,  thus  leading  to  a 
different  divergence  factor  for  each  of  the  two  different  characteristic 
ray  paths. 


The  back- scattering  geometrical  divergence  factor  for  the  "v"  ray  path 
is  given  by 


(6-i°>  o„r 

which  can  be  simplified  to  give 

(6.11)  DiRv  =  * 

T-v-K] 


XtAn.  I  ~  CO&.  "2.  A  &L\r 

Afii ~o 


The  back- scattering  geometrical  divergence  factor  for  the  Mvds”  ray 
path  is  given  by 


(6-12)  .  -t 

_  ?•  °H  JEi/m.  I  *  A  A 

tgiTo  i  iv 


(■ 

>! 

li 


A  0iyr  +  m.  okc  ccL  €v+Aju»v  A5iw.  +  ^  +  ( €f- Aim*  A  Q& 


^1-  C*4  2. 

fcl'"' 
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'4  ii  4 


2  ( Aft,)1]  *j 
|er+WAft, +[(ft*  l^-^Aft,)1]^]-^]  ^Aft^rJ  yi||  | 


which,  thru  the  application  of  l'Hospital's  rule,  can  be  evaluated  to 
give 


(6. 13) 


4  R<v, 

Mi1 


I  +  na  —  , — 

>  [e. 


VT 


♦  er\.<rr*l  *]K 

J  > 


Using  the  relative  phase  factor  /3r  ,  as  defined  in  Appendix  1 , 

simplifies  equation  (6. 13)  to 


(6. 14) 


Di.  Rvats  n  ~ 


4RR 


I  +  /wt 


z 


(1)  S.  Silver,  Microwave  Antenna  Theory  and  Design,  McGraw-Hill, 
1949,  pp.  HT-  114.  " 
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APPENDIX  A-7 

GEOMETRICAL  OPTICS  ELECTROMAGNETIC  WAVE 
BACK-SCATTERING  CROSS  SECTION  OF  A  UNIFORM 
LOSSY  DIELECTRIC-COATED  METALLIC  SPHERE 


Detailed  here  is  the  theoretical  determination  of  the  electromagnetic 
wave  back- scattering  cross  section  (radar  cross  section)  of  a  per¬ 
fectly  reflective  sphere  coated  with  a  uniform  layer  of  lossy  dielectric 
material  via  the  geometrical  optics  approach  using  the  ray  path  field 
intensity  addition  technique. 


Using  the  ray  path  field  intensity  addition  technique,  the  general 
expression  for  the  back- scattering  cross  section  (from  Appendix  3) 
is  given  by 


(7.1) 


a  -  ak  r, 


TT  /%r 


v*i 


T* 

TT  T« 


t«t 


+•1 


I 


e 


z 


The  sum  over  k  can  be  expanded  into  two  sums,  each  over  one  of  the 
two  different  characteristic  ray  paths  (”vM  and  "vdsM;  see  Appendix  5). 
This  is  reasonably  accurate  when  the  coating  thickness  (h)  is  small 
compared  to  the  sphere  radius  (r)  because  the  effects  of  encircling  ray 
paths  are  negligible  in  this  case. 
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Inspection  of  the  characteristic  ray  paths  shows  that  the  sum  of  all  the 
contributions  from  the  "v"  characteristic  ray  path  reduces  to  a  single 
term,  while  the  sum  of  the  contributions  from  the  "vds"  characteristic 
ray  path  consists  of  an  infinite  number  of  terms.  Expanding  the  sum 
over  k  into  sums  over  the  two  different  characteristic  ray  path 
contributions  gives 


(7.2)  <T=  4  IT  Rr 


-'/a 

/°vrd.  DlRv 


“i  W  \T 


+  I 

mi*  l 


D 


-54 

iRvrcLsrm 


rr  V*  TT  TV^t 

tz\ 


Bvr^Snm.  -j/3^5^  L 4, 

TT  e 
1*» 


Avglsrm  ^-vrdliwia 

TT  c 


*.»  » 


Examination  of  the  "vds”  characteristic  ray  path  shows  that 


(7.3) 

TT  =  (-lP 

r  =  i 

r»rv  —  1 

t/°^) 

(7.4) 

TvrdLtmv 

TJ  =  Trl 

(7. 5) 

TT  c  = 

1 

C  c 
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A  —  Of  v-dLso.  L  "  ^flC  D 

(7.6)  XT  =  e 

CL  *  I 

Appendix  4  shows  for  normal  incidence  that 


(7-7)  /°JLr  =  "/°vri 
(7.8)  'HriTi,  = 

and  Appendix  5  shows  for  the  back-scattering  condition  that 


(7.9) 

(7.10) 

(7.11) 


L-Vnt  =  ° 
VvrdLSfm 
^’Ov'JlSmv 


o 

2.  onrv  A. 


Using  Equations  (7. 3)  through  (7. 11)  in  equation  (7. 2)  gives 


i 

(7.12)  <f  =  4TT  R| 


/°vd.  D 


-/«. 

iR\r 


_  ji^Vaj  r 

^  ^iRvr«L*nm, 

< \rd-  mv  =  l 


^4  <=  e  j 


41V 
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The  explicit  forms  of  the  various  factors  in  equation  (7. 12)  are 
necessary  for  completion  of  the  formula. 

Appendix  A-4  shows  for  normal  incidence  that 


(7.13)  /^-L 


•- 1 

-  i  ] 

i  /z 

1  +  1 

-j«r  ) 

and  Appendix  A-6  shows  for  the  back- scattering  condition  that 


(7. 14) 

-!4 

DLRvr  " 

1  1  1  1 

(7.15) 

-  ‘/at 

DiRgrASnwv 

_  (n-ft.)  j 

**  *R 

f.  .  All 

\  -V-  rm. - — — 

*•  /®r 

-1 

It  can  be  shown  ^  that  the  phase  factor  (  /&&  )  and  attenuation  factor 
(CCJl  )  for  a  plane  electromagnetic  wave  propagating  normally  into  a 
lossy  dielectric  medium  are  the  same  as  those  derived  for  plane 
electromagnetic  wave  propagation  in  the  unbounded  lossy  dielectric 
medium.  These  factors  are 


(7. 16)  /Sj 
<7.!7) 


V5~  TT 

Xvr 

VF  TT 

\>r 


[e,  ♦  (C  +  < ) 
[_e-  ♦  {«  *  -v  < 


'/z 
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The  relationship  between  the  electromagnetic  wave  wavelength  in  the 
lossy  dielectric  medium  (  X  ^  )  and  the  free-space  wavelength 
(  \  v  )  is 


(7.18) 


2.TT 


&  \r 


'/z 


Vz  ~ 


Using  equations  (7. 14)  and  (7. 15)  in  equation  (7. 12),  and  normalizing 
with  respect  to  the  radar  cross  section  of  the  uncoated  sphere 
(  IT  r 2  ),  gives 


(7.19)  — 2— 

TT  rz 


( 1  —  sd ) 
(°  wd. 


e 


Z 


where  vrJL  >  ,  and  are  given  by  equations  (7.13), 

(7. 16),  and  (7. 17)  respectively. 


Note  that  when 


/Vjl  - 

(7.20) 


oo 

l 

nnel 


- 

£  1 

-i 

L  r 

»  • 
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the  radar  Gross  section  of  the  uniform  lossy  dielectric-coated 

) 

metallic  sphere  will  vanish. 

A  useful  approximation  to  equation  (7. 19)  for  the  radar  cross  section 
of  a  uniform  lossy  dielectric-coated  metallic  sphere  can  be  made  if 
the  condition 


(7.21)  <<  I 

**  /Of 


is  met.  Under  this  condition,  the  approximation 


(7.  22) 


I  + 


A  1 

* 


rm. 


can  be  made,  which  modifies  equation  (7. 19)  to  the  form 


(7.23) 

<r  ~  i 

ITr*  ~  l 

-4 r 

oo 

/°vi  «  <=• 

-OCjL  2-A 

L 

wn.-\ 
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Using  the  formula 


(7.24)  Y  X 
mv«i 

the  sum  in  equation 


n  _  X 

I  -  X 

(7.  23)  can  be  closed,  giving 


(7.25) 


*  1  + 


/°vrd.  * 


-j/Sj.ail  -ocA 

e. 


\ 


-ft.  i 

\r 


I 


-tit*. 


t 


*  -#•  1 


Note  that  when 


(7.26) 


ircL 


-  Kj.  Z  A. 


I  + 


A  I 
►  £»» 
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the  approximate  radar  cross  section  of  the  uniform  lossy  dielectric- 
coated  metallic  sphere  will  vanish.  Since  this  is  a  conditional  equation 
in  complex  variables  it  represents  two  equations  in  real  variables 
defining  the  necessary  relationships  between  the  dielectric  constant 
(  £r  ),  the  relative  conductivity  (  G  y.  ),  the  coating  thickness 
(  ^  ),  the  sphere  radius  (  T  ),  and  the  electromagnetic  wave  free- 
space  wavelength  (  X  ^  )  that  must  be  satisfied  in  order  for  the  radar 
cross  section  as  given  by  equation  (7.  25)  to  vanish.  Equation  (7.  26) 
can  be  separated  into  the  two  conditional  equations 


(7.27) 


TT  4-  oxc.tam. 


4TT  t,  „  2 
i  x  ”  OCjL  ~  ftcL 
_ 

4TT  cC^ 


where  k  =  1,2,3 . and  indicates  the  multiplicity  of  values  for  which 

the  radar  cross  section  can  vanish.  Equation  (7.27)  gives  the  necessary 
relationship  between  the  electromagnetic  properties  of  the  dielectric 
material  and  equation  (7.  28)  gives  the  necessary  coating  thickness. 
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The  formulas  for  the  detailed  variation  of  the  radar  cross  section  of  a 
uniform  lossy  dielectric-coated  metallic  sphere  (equation  (7.19)  and 
equation  (7.25)  for  the  exact  and  approximate  formulations,  respec¬ 
tively)  are  sufficiently  complex  to  preclude  understanding  of  the 
detailed  functional  variations  involved  merely  by  inspection  of  the 
equations.  It  is  therefore  necessary  to  numerically  calculate  and 
plot  the  radar  cross  section  in  a  suitably  normalized  form  in  order 
to  gain  insight  into  the  detailed  functional  behavior. 

Using  equation  (7. 19),  the  radar  cross  section  of  a  uniform  lossy 
dielectric-coated  metallic  sphere  has  been  calculated  for  certain 
combinations  of  the  various  parameters  involved.  The  results  of  these 
calculations  are  shown  in  Figures  A7-1  and  A7-2,  where  the  norma¬ 
lized  radar  cross  section  is  plotted  as  a  function  of  the  dielectric 
constant  with  the  normalized  coating  thickness,  normalized  wave¬ 
length,  and  relative  conductivity  as  parameters.  The  normalized 
coating  thickness  (  «-A  )  has  the  same  value  (  o.l  )  for  all  the 
curves  shown  in  Figures  A7-1  and  A7-2.  The  normalized  wavelength 
(  X*  / r  )  values  (0. 034215  and  0. 25146  for  Figures  A7-1  and 
A7-2,  respectively)  are  chosen  to  provide  the  longest  wavelengths 
(k  =  1)  that  cause  zero  radar  cross  section  to  occur  for  zero  relative 
conductivity  (in  Figure  A7-1)  and  for  a  non-zero  (0.62772)  relative 
conductivity  (in  Figure  A7-2),  respectively.  Several  values  of  rela¬ 
tive  conductivity  are  used  to  show  the  transitional  behavior  from  the 
lossless  to  the  highly  lossy  condition. 

The  geometrical  optics  approach  to  the  determination  of  the  radar 
cross  section  of  a  uniform  lossy  dielectric-coated  metallic  sphere  is 
applicable  and  valid  when  the  radius  of  the  sphere  is  much  larger  than 
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47-11 


Figure  A7-2  Radar  Cross  Section  of  a  Uniform  Lossy  Dielectric- 
Coated  Metallic  Sphere,  X  /r  =  0.  25146 
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the  wavelength  of  the  incident  electromagnetic  wave  in  free  space  and 
in  the  dielectric.  Using  equation  (7. 18)  these  applicability  criteria 
become 

(7.29)  -  <<  r 

and 

(7.30)  -4*-  «  V 

Pr 


These  geometrical  requirements  limit  the  applicability  of  the  geo 
metrical  optics  approximation  when  the  dielectric  constant  and  'he 
relative  conductivity  both  closely  approach  zero. 


(1)  J.  A.  Stratton,  Electromagnetic  Theory,  McGraw-Hill,  1941, 
pp.  500-524. 
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4  APPENDIX  A-8 

GEOMETRICAL  OPTICS  RADAR  CROSS  SECTION  OF  A 
METALLIC  SPHERE  COATED  WITH  A 
UNIFORM  LOSSY  PLASMA 

The  radar  cross  section  <T  ,  obtained  by  the  geometrical  optics 
approach,  of  a  perfectly  reflective  (metallic)  sphere  of  radius  Y* 
covered  with  a  uniform  coating  of  lossy  plasma  is  given  here.  The 
plasma  is  assumed  to  have  uniform  thickness  all  around  the 
sphere  and  no  spatial  variation  of  its  electromagnetic  properties 
throughout  the  coating.  Furthermore,  it  is  assumed  that  the  plasma 
can  be  characterized  (electromagnetic ally)  by  a  permeability 
a  capacitivity  ,  and  a  conductivity  analogous  to  the 

situation  for  non-ionized  media.  Thus,  the  situation  becomes  that  of  a 
uniform  lossy  dielectric-coated  metallic  sphere.  The  rigorously 
obtained  expression  for  the  radar  cross  section  of  this  composite 
structure  is  given  in  Appendix  2,  with  appropriate  substitutions  for 
the  dielectric  properties  of  the  plasma  from  Appendix  1.  However, 
some  of  the  general  characteristics  of  the  behavior  of  the  radar  cross 
section  under  various  parameter  variations  can  be  better  seen  from 
the  geometrical  optics  approach.  Using  the  results  from  Appendixes 
1  and  7,  the  normalized  radar  cross  section  is  given  by 
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Jj 

where  D  is  the  relative  geometrical  divergence  factor  for  the  nm.  ray  path 
m 

thru  the  plasma  coating,  Hp,  is  the  reflection  coefficient  for  a  plane  transverse 
electromagnetic  wave  propagating  normally  from  free  space  into  a  plane  inter¬ 
face  between  free  space  and  a  uniform  plasma  having  the  electromagnetic  pro¬ 
perties  of  the  coating,  is  the  propagation  factor  for  a  plane  transverse 

electromagnetic  wave  propagating  thru  an  unbounded  uniform  plasma  having  the 
electromagnetic  properties  of  the  coating,  and  fir  is  the  relative  phase  factor 
for  the  same  situation.  Appendixes  1,  4  and  7  give 


(a  2) 


(8.3) 


(8.4) 


(8.5) 


1  4  mx 


A 


-1 


A 


rj 

[*«*♦ 

yM 

’/»  o  *  . 

♦  + 

i  *■ 

1  -  -ILCX 

where 

* 

*4  * 

CO  = 


cUp  / co  =  normalized  plasma  frequency, 

Vc  /<*>  =  normalized  electron  collision  frequency, 

angular  plasma  frequency  (  2.  IT  ), 
electron  collision  frequency, 
angular  radar  frequency  (  2.  TT  -f  ),  and 
free-space  wavelength  of  the  radar  wave. 
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It  can  be  seen  from  equation  (8. 1)  that  when  the  plasma  coating  is 
highly  reflective  (  £  i  )  the  radar  cross  section  4  approaches 

the  value  TTr*(i  + &/»•)*,  which  is  the  value  for  a  perfectly  reflective 
sphere  having  a  radius  equal  to  the  outer  radius  of  the  plasma  coating. 
When  the  plasma  coating  is  negligibly  reflective  (  J  -  O  )  the 
radar  cross  section  <T  approaches  the  value  TT  r16,  ,  which  is  the 
value  for  the  bare  metallic  sphere,  as  would  be  expected. 


Note  that  when 


(8.6) 


n.  = 


2\  00 

-jA  I  D 

n  Unrs. 


f> 


rm,=  t 


rv  e 


-iv  \  n"v 


the  radar  cross  section  will  vanish.  This  condition  is  possible  but  is 
very  difficult  to  express  concisely  in  terms  of  the  variables  Y  ,  , 

X  v  ,  £1  ,  and  Xle  because  of  the  infinite  sum  involved. 


A  useful  approximation  to  equation  (8. 1)  can  be  found  (see  Appendix  7) 
when 


(8.7) 


«  \ 


Under  this  condition  the  infinite  sum  in  equation  (8. 1)  can  be  closed 
approximately,  giving 

p  -42^  2 

r>  -  d  e  r 

-Hu**. 

i-rfDer 

where 


(8.8) 


TTr* 


4  1 
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(8.9) 


D 


Equation  (8.8)  has  the  same  behavior  as  equation  (8. 1)  for  both  highly 
reflective  and  negligibly  reflective  plasma  coatings.  Also,  when 


(8. 10) 


-KZ*. 

D  e  * 


the  radar  cross  section  will  vanish.  This  can  occur  only  when  both  the 
magnitudes  and  the  phase  angles  of  the  two  terms  are  equal.  These 
two  conditions  give  the  parametric  requirements  for  vanishing  of  the 
radar  cross  section,  and  can  be  written  in  the  form 


fir  0_  I"  **•*  +  1 1  +  y®f)  1  -.a. 

**  +  JI-A.)  *  J 


(8. 12) 


A  - 

\\r 


Wr 


4  a*£.t*/vv 


I  -  OOr  - 


2X, 


where  ■ft.  =  1,2,3,....  and  indicates  the  multiplicity  of  values  for  which 
the  radar  cross  section  can  vanish,  and  <C  y>  is  the  relative  attenuation 
factor  given  by 


(8.13) 
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Equation  (8. 11)  gives,  indirectly,  the  necessary  relationship  between 
the  electromagnetic  properties  of  the  plasma  and  equation  (8. 12)  gives 
the  necessary  plasma  thickness  required  for  the  radar  cross  section 
to  vanish.  In  a  three-dimensional  space  whose  coordinates  are  , 
Sic  ,  and  Jk/xr  the  conditions  for  vanishing  radar  cross  section 
are  lines,  one  line  for  each  value  of  &  .  It  is  necessary  to  solve 
equations  (8. 11)  and  (8. 12)  numerically  in  order  to  see  the  behavior  of 
these  conditions  in  terms  of  the  pertinent  variables.  Figures  A8-1 
and  A8-2  show  a  typical  solution.  Note  that  the  conditions  given  by 
both  figures  must  be  satisfied  simultaneously.  It  can  be  seen  from 
these  figures  that  the  normalized  plasma  frequency  (cOf»/co  )  re¬ 
quired  for  zero  reflection  is  of  the  order  of  unity  for  small  (less  than 
unity)  values  of  normalized  electron  collision  frequency  (tJc/co)  and 
varies  approximately  as  the  square  root  of  the  normalized  electron 
collision  frequency  for  large  (greater  than  unity)  values  of  normalized 
electron  collision  frequency.  It  is  important  to  note  that  the  minimum 
plasma  layer  thickness  (A.  )  for  which  zero  reflection  can  occur  is  of 
the  order  of  one-third  of  the  free- space  wavelength  of  the  radar  wave. 
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Figure  A8-2  Conditions  for  Vanishing  Radar  Cross  Section 
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APPENDIX  A-9 

ELECTROMAGNETIC  WAVE  REFLECTION  AT  NORMAL 
INCIDENCE  FROM  AN  INFINITE  METALLIC  FLAT  PLATE 
COATED  WITH  A  UNIFORM  LOSSY  PLASMA 


Derived  here  is  the  reflection  of  a  plane  transverse  electromagnetic 
wave  normally  incident  on  an  infinite  perfectly  reflective  (metallic) 
flat  plate  coated  with  a  uniform  lossy  plasma.  The  plasma  coating  is 
assumed  to  have  a  uniform  thickness  h,  and  no  spatial  variation  of 
its  electromagnetic  properties  throughout  the  layer.  Furthermore,  it 
is  assumed  that  the  plasma  can  be  completely  characterized  electro- 
magnetically  by  an  equivalent  permeability,  dielectric  constant,  and 
conductivity  as  given  in  Appendix  A-l. 

In  the  free-space  region  outside  the  plasma  layer  the  incident  electro¬ 
magnetic  wave  is  given  by  the  electric  field  (  )  distribution 

__  .  _  j  cot  -Vy  X 

(9. 1)  »  1*  Ei  ^ 

and  the  associated  magnetic  field  (  {Hi,  )  distribution 

j  /I —  j  o»  t  -  r 

(9.2)  .  IHi.  -  Ei  e 

where  the  wave  is  taken  to  be  propagating  in  the  positive  z  direction 
(normal  to  the  flat  plate)  and  has  the  electric  field  vector  polarized  in 
the  positive  x  direction.  The  propagation  factor  in  free  space  (  X  ^  ) 
is  given  by 
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(9.3)  tr  -  j  ^r~ 

A  y- 


where 


1 


%  unit  vectors  in  the  x  and  y  directions 

E  L  =  amplitude  of  the  electric  field  of  the  incident 

CO  =  angular  frequency  of  the  incident  wave  (  2.V( 

t  =  time 

=  free-space  wavelength  of  the  incident  wave 
€y-  =  capacitivity  of  free  space 
Jl\r  =  permeability  of  free  space. 


wave 

f  ) 


In  the  free-space  region  outside  the  plasma  layer  the  reflected  wave 
is  given  by  the  electric  field  (  {£  t  )  distribution 


(9.4) 


j  <■>*  +<V 

E,  =  l*Er<= 


and  the  associated  magnetic  field  (  (H  r  )  distribution 


(9.5) 


In  the  plasma  layer  the  total  field  consists  of  forward-and  backward¬ 
traveling  waves.  The  forward-traveling  wave  is  given  by  the  electric 
field  (  CE  ^  )  distribution 
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i  jwt  -  /fZ. 

(Ef  ”  lx  €2 

and  the  associated  magnetic  field  (  IH*  )  distribution 
(9'7) 

The  propagation  factor  in  plasma  (  tflja.  )  is  given  by 

(9-8)  t¥  «  j^7K-jtfV)'/t 

where 

£r  =  dielectric  constant  of  plasma 
OV  =  relative  conductivity  of  plasma 

The  backward-traveling  wave  is  given  by  the  electric  field  (  (E  ) 
distribution 

jcot  f  Z- 

(9.9)  E*  =  1*  E±  e 

and  the  associated  magnetic  field  (  )  distribution 

The  reflection  coefficient  (  P  )  for  the  composite  structure  is 
defined  as 
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at  the  plasma  surface 

Assuming  the  plasma  surface  to  be  at  z  =  0  and  the  metallic  flat 
plate  to  be  at  z  =  h  gives 

-  Er 

(9.12)  r  = 

Ei. 

Applying  the  electromagnetic  boundary  conditions  at  the  plasma  surface 
of  continuous  tangential  electric  field  gives 

(9.13)  Et  +  =  Ef  ♦ 

and  continuous  tangential  magnetic  field  gives 

(9.14)  Ei  -Et  =  (6,-j«;)!4|Ef-Ef| 

Applying  the  electromagnetic  boundary  conditions  at  the  metallic  flat 
plate  surface  of  zero  tangential  electric  field  gives 


Solving  equations  (9. 13)  thru  (9. 15)  gives 
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i  -  (£,  -  Jj;)14  _ 

'  *  K--K)  1/4  ~ 

-3V«-R.  ,.Uv-i<rr) ’/*■ 

i+k-Kl* 


and  using  the  definition  of  the  reflection  coefficient  for  a  plane  inter¬ 
face  between  free  space  and  an  unbounded  plasma  (  P^,  )  given  from 
Appendix  4  by 


(9. 17) 


1  - 

-  j  <r^ 

lA 

1  ♦  i 

^  r  "  j 

Vi 

gives  the  overall  reflection  coefficient  (  P  )  as 

ri  fjt,  **  ® 

(9-i8)  r  =•  — * - j?rrr 

i  -  r>  e  f 

It  can  be  seen  from  equation  (9. 18)  that  when 


(9. 19) 
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there  will  be  no  reflected  wave.  Using  the  results  of  Appendix  1  which 


give 

(9.  20) 

P  — 

I  - 

,  .A 

1  1  l-jilc  1 

^  ’ 

1  + 

1  -  J>v)  A 

and 

(9.21) 

^  - 

1  ttr 

1  / 

where 

•Q-p  = 

fic  =\Jc/co 
=  2^ 

>4  = 

oj  =2irf 


normalized  plasma  frequency 
normalized  electron  collision  frequency 
angular  plasma  frequency 
electron  collision  frequency 
angular  frequency  of  incident  wave 


the  plasma  properties  necessary  for  vanishing  reflection  can  be 
determined  from  equation  (9. 19).  Since  equation  (9. 19)  involves 
complex  variables  it  is  equivalent  to  two  ordinary  algebraic  equations 
which  can  be  written  in  the  forms 


(9. 22) 


(t*- )  - 


JL. 

-1 

ex,, 

-  *  * + 

I'-A 
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and 

(9.23) 

A  _ 

1 

X\r 

4T Tflr 

where 

k  - 

i,  2. 

(9.  24) 

<r,=  j 

(9.  25) 

A  4- 

[-Ac* 

] 


x  «  2 


and 


2  (l  +  n,?) 

!4 


Z  1'-^]  J 


Equation  (9.22)  gives,  indirectly,  the  necessary  relationship  between 
Xl^  and  XI  c  ,  and  Equation  (9.  23)  gives  the  necessary  relation¬ 
ship  between  A,  and  Xy  for  vanishing  reflection.  The  multiple 
values  of  &  indicate  the  multiplicity  of  thicknesses  for  which  the 
vanishing  condition  can  be  obtained.  In  a  three-dimensional  space 
whose  coordinates  are  •Qf ,  -n-c  ,  and  a/xr  the  conditions  for 
zero  reflection  are  curved  lines,  one  line  for  each  value  of  k.  It  is 
necessary  to  solve  equations  (9.  22)  and  (9.  23)  numerically  in  order  to 
see  the  behavior  of  these  conditions  in  terms  of  the  pertinent  variables. 
Figures  A9-1  and  A9-2  show  the  solution.  Note  that  the  conditions 
given  by  both  figures  must  be  satisfied  simultaneously  for  some  given 
k  value.  It  can  be  seen  from  these  figures  that  the  normalized  plasma 
frequency  (  -a*  )  required  for  zero  reflection  is  of  the  order  of  unity 
for  small  (less  than  unity)  values  of  normalized  electron  collision 
frequency  (  XI  c  )  and  varies  approximately  as  the  square  root  of  the 
normalized-electron  collision  frequency  for  large  (greater  than  unity) 
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Figure  A9-2  Conditions  for  Zero  Reflection  ( 
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values  of  normalized  electron  collision  frequency.  It  is  important  to 
note  that  the  minimum  plasma  layer  thickness  (  )  for  which  zero 

reflection  can  occur  is  of  the  order  of  one-third  of  the  free-space 
wavelength  of  the  incident  electromagnetic  wave. 
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APPENDIX  A-10 

RAY  PATHS  IN  NON-UNIFORM  LOSSY 
DIELECTRIC  MEDIA 


The  concepts  of  geometrical  optics  (ray  theory)  are  often  very  useful 
in  the  general  consideration  of  the  propagation  of  electromagnetic 
waves,  and  sometimes  can  be  applied  to  greatly  simplify  the  method 
of  analytical  solution  of  electromagnetic  wave  propagation  problems, 
when  they  are  applicable.  In  general,  the  concepts  of  geometrical 
optics  are  applicable,  and  the  results  of  solutions  using  these 
concepts  are  reasonably  accurate,  when 

(1)  the  fractional  change  in  the  constitutive  parameters  of  the 
media  thru  which  propagation  occurs  is  small  over  a  distance 
of  one  wavelength  in  the  direction  of  propagation  (for  con¬ 
tinuously  varying  constitutive  parameters)  or  if  the  total 
change  in  tne  constitutive  parameters  takes  place  over  a 
distance  very  much  smaller  than  one  wavelength  (for  an 
abrupt  boundary  between  two  regions),  and 

(2)  the  fractional  change  in  the  spacing  between  adjacent  rays 
over  a  distance  along  the  rays  of  one  wavelength  is  small.  This 
means  that  the  geometrical  optics  solution  is  open  to  question 
if 

(a)  the  rays  pass  thru  a  focus,  or 

(b)  the  rays  are  reflected  or  refracted  at  a  surface  with  a 
size  or  radius  of  curvature  that  is  not  large  compared  to 
the  wavelength  . 
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A  ray  path  can  be  considered  to  be  the  path  followed  by  a  narrow 
pencil  beam  of  electromagnetic  radiation  as  it  propagates  thru  a 
medium.  In  uniform  media  ray  paths  are  straight  lines.  In  non- 
uniform  media  ray  paths  are  curved  toward  the  direction  of  decreasing 
velocity  of  propagation.  The  general  vector  differential  equation  that 
defines  a  ray  path  in  a  non-uniform  isotropic  lossless  medium  is  ^ 


(10.1) 


where 

S  =  distance  along  the  ray  path  from  some  reference  point  on 
the  ray  path 

O'*  =  position  vector  from  some  reference  origin  to  the  ray  path 

m.  =  index  of  refraction  of  the  medium 
Y  =  gradient  operator. 


In  a  non-uniform  isotropic  lossy  medium  the  index  of  refraction  as 
usually  defined  is  complex.  In  this  case  the  index  of  refraction  is  re¬ 
placed  by  the  electromagnetic  wave  relative  phase  factor  ,  which 
is 


(10.2)  fif  - 

where 

**\r  =  wavelength  of  the  electromagnetic  wave  in  free- space 
A.  =  wavelength  of  the  electromagnetic  wave  in  the  medium 
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thus  giving  the  ray  path  equation  in  a  non-uniform  isotropic  lossy 
medium  as 


In  a  lossy  dielectric  material  with  dielectric  constant  €y»  and  con¬ 
ductivity  ,  the  relative  phase  factor  is 

(10. 4)  /9,  --  ~  [ef  +  I  €*  ♦  <rr* ) 14  ] /l 

where 

C\ i  =  £  -gj  =  relative  conductivity  of  the  dielectric, 

Gy.  =  capacitivity  of  free  space,  and 

=  ZtTf  =  angular  frequency  of  the  electromagnetic  wave. 

It  should  be  noted  that  the  wavelength  of  an  electromagnetic  wave  in  a 
lossy  dielectric  material  depends  on  the  conductivity  as  well  as  on  the 
dielectric  constant  of  the  dielectric  material. 

Two  specific  geometrical  configurations  of  interest  here  are  the 
circular  cylindrical  and  the  spherical  coordinate  systems. 

Consider  first  the  circular  cylindrical  coordinate  system,  shown  in 
Figure  A10-1.  Assume  that  in  this  system  /3r  may  be  a  function  of 
both  /°  and  z.  but  not  of  ^  .  Now  ray  paths  entering  the 

region  traveling  in  a  =  constant  plane  will  remain  in  that  plane, 

although  they  may  bend  toward  or  away  from  the  z.  -axis  as  they 
progress  in  the  Z  direction.  In  this  case,  the  position  vector  from 
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the  origin  to  a  point  on  the  ray  path  (“•ff*  )  is  *-■ 

(10.5)  p  =  \fi/o  ♦  1*2. 

where  1  p  and  i  *  are  unit  vectors  in  the  /°  and  z.  directions, 
respectively.  An  incremental  length  along  the  ray  path  (  is  )  is 
given  by 

(10.6)  (is)  *  =  (i/0)*  ♦  U*)* 


and  the  gradient  of  the  relative  phase  factor  (  V/9r  )  is  given  by 


(10.7) 


j  ,  A 

lfi  Tp  lz-  IT 


In  order  to  describe  the  ray  path  it  is  necessary  to  write  equation 
(10.3)  in  a  form  such  that  it  expresses  the  relationship  between  /© 
and  Z.  along  the  ray  path,  independent  of  the  arc  length  (s).  The 
chain  rule  of  differentiation  shows  that 

(10.8)  _i_  -  ^ 

is  “  is  i/° 


and  equation  (10.  6)  gives 


TR63-217A 


Using  equations  (10.8)  and  (10.9)  gives 


(10.10) 


dl T  - 
els 


I  + 


cl  Z 
2.  d/O 


Thus  equation  (10. 3)  becomes 

,-34 


^/®r 

dz 


Separating  this  vector  equation  into  its  two  scalar  equations  gives 


(10.12) 


and 


(10.13) 


Expanding  the  derivative  in  equation  (10. 13)  and  substituting  equation 
(10. 12)  into  it  gives 


(10.14) 


cU  J_ 

a/0  P, 


I  cLx]Z  » 

l  JL/°  I  3z- 


* 


J _ 

ft  d/° 
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In  order  to  find  the  ray  path  it  is  necessary  to  solve  equation  (10. 14) 
for  z.  in  terms  of  P  . 


Consider  the  special  case  where  /3r  is  not  a  function  of  zl  .  Then 
equation  (10. 14)  simplifies  to 


do.  is)  + 1 Ah.  |3 -L-  -  0 

JL/>*  dL/®  /8,  U/®)  fir  d/° 


This  is  a  Bernoulli  equation  for 
give 


dlz. 

±P 


which  can  be  solved 


(2) 


to 


(10.16) 


d/° 


where  K  is  a  constant  of  integration  and  the  algebraic  sign  determines 
the  direction  of  propagation  along  the  ray  path.  Equation  (10. 16)  can  be 
integrated  to  give 

(10.17)  Z-Z0  =  ±  j  VXd.f> 

1 

where  /£  ;  Z0  is  some  initial  point  on  the  ray  path  and 

is  a  function  of  p  .  Note  that  the  constant  of  integration  (K)  can 
be  expressed  as 

(10.18)  K  =  ^  1  + 
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Thus,  in  order  to  define  a  particular  ray  path  it  is  necessary  to 
specify  the  location  of  one  point  on  the  ray  path  and  the  slope  of  the 
ray  path  at  that  point.  Equation  (10. 18)  can  be  put  in  the  form 

(10. 19)  K  =  -3TT  (  I  *  K  ) 

'o 

where  ^  is  the  angle  between  the  z.-axis  and  the  tangent  to  the 
ray  path  at  the  initial  point  /£  ;  Z0  . 

Consider  the  particular  form  of  variation  of  /Qr  given  by 

(10.20)  |  +  I/O1) 


which  is  parabolic  in  /&tz  .  When  is  algebraically  positive 
then  /Gy  increases  as  /°  increases.  In  this  case  the  ray  paths 
are  bent  away  from  the  z.  -axis  (de-focused).  When  is  alge¬ 
braically  negative  then  /6r  decreases  as  increases.  In  this 
case  the  ray  paths  are  bent  toward  the  z.  -axis  (focused).  Substitution 
of  equation  (10.20)  into  equation  (10. 17)  and  integration  gives 


,  .  z-z 

(10.21) 


when  ir  is  algebraically  positive,  and 

z~z0  Z/9','  l_*K)  /a  “*•**'"•  (  [  I  -  K-I  /3 ] 
(10.22)  1  5 

[/» ( -/2‘)  -/S  ( -/•*  ] 
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when  fir  is  algebraically  negative.  The  general  behavior  of  these  two 
different  types  of  ray  paths  can  be  demonstrated  by  considering  the 
special  case  where  is  zero.  Then  equations  (10.  21)  and  (10.22) 

reduce  to 


(10.23)  Z--Zo-0~s  '  \W)/l 


and 


(10.24)  Z-ZQ  = 


=  A,  Mk)  '  OKC-Xi/n.  |/°(|.K_1/a^i.)  J 


respectively.  These  can  be  rewritten  in  the  form 

—i  a 


(10.25) 


fi  -  (frx)'4  |z-z.)] 


and 


(10.26)  /°  =  (-  |'4K)^(Z-Z.)] 


which  show  that  as  the  ray  travels  in  the  Z.  direction  /°  increases 
monotonically  (de-focusing)  when  -fir  is  algebraically  positive  and 
increases  and  decreases  periodically  (periodic  focusing)  when  &  is 
algebraically  negative. 


Now  consider  the  spherical  coordinate  system,  shown  in  Figure  A10-2. 
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Assume  that  in  this  system  may  be  a  function  of  both  R  and 

Q  but  not  of  $  .  Now  ray  paths  entering  the  region  traveling  in 

a  =  constant  plane  will  remain  in  that  plane,  although  they  may 

bend  toward  or  away  from  the  origin  as  they  progress  thru  the  region. 
In  this  case,  the  position  vector  from  the  origin  to  a  point  on  the  ray 
path  (  O'*  )  is 

(10.27)  0*  =  1R  * 

where  is  a  unit  vector  in  the  R  direction.  An  incremental  length 
along  the  ray  path  (  els  )  is  given  by 

(10.28)  d.S  Z  =  (<±R)X  +  ft*  (<±0)Z 

and  the  gradient  of  the  relative  phase  factor  (  VA  )  is  given  by 

(10-29)  VA  =  i«-f£  * 

where  f  ^  is  a  unit  vector  in  the  Q  direction. 

In  order  to  describe  the  ray  path  it  is  necessary  to  write  equation 
(10. 3)  in  a  form  such  that  it  expresses  the  relationship  between  R.  and 
Q  along  the  ray  path,  independent  of  the  arc  length (s).  The  chain 
rule  of  differentiation  shows  that 


(10.30) 


d  s  _aL_ 

”  el  5  JiQ 


and  equation  (10. 28)  gives 


(10.31) 
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Using  equations  (10. 30)  and  (10. 31)  gives 


(10.32) 


dor 

els 


-  >**  i#r]  i* 


R  iff 


Thus  equation  (10.3)  becomes 


(10. 33) 


R*+i£  4 
R  +  60 


j  i 

1a  r 


==  iri&* 


R  a*.  *  a# 

Separating  this  vector  equation  into  its  two  scalar  equations  gives 


(10. 34) 


M#rr  i*i^i#ri 


Ufl*  I 


60  60 


M 


a.  (oLR 

R  +  iff, 


’ll  1  ■ 


dte 


and 


Equations  (10. 34)  and  (10. 35)  can  be  combined  to  give 


(10. 36) 


dzR 

<LR  i 

dfif  |<±R)* 

Z  1  $  fir  \ 

el  0* 

6.8  fit 

d0  Ufll 

R  /S^T  3R  1 

Ji£l3J _ L  ±0*  -  r  .  r>*  i  »£>• 

*U0  ifi  fir  dR. 
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In  order  to  find  the  ray  path  it  is  necessary  to  solve  equation  (10.36) 
for  R  in  terms  of  6  . 


Consider  the  special  case  where  the  region  of  interest  is  confined  to 
the  spherical  shell  between  two  closely  spaced  concentric  spherical 
boundaries  of  radius  and  R^.  (where  RX>R,,  and 
and  where  within  this  region  is  not  a  function  of  R  .  Equation 
(10. 36)  reduces  to 


(10.37) 


dR  *  I  olRi*  Z 

A6Z  AB  0*  <L6  l  del  R 

.  I A  R  \3  i  J _ A  Py  o 

*R*  fiy  dL 9 


0 


Transforming  to  new  variables  p  and  Z.  ,  defined  by 


(10.38)  (O  =  R*  8 
and 

(10.39)  z.  -  R»-  R 


changes  equation  (10. 37)  to  the  form 

(10.40)  A*Z.  AZ  J_  AjPt  Uz]8,  2. 

4.1°*  Ap  fir  Ap  l  Apt  R t-Z. 


Ux|s  i 

A0f 

U p]  Pr 

d/° 
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and  when  is  very  large  and  z.  remains  small  this  reduces  to 


(10. 41) 


dZ  I  dfir  + 

d/°  /6r'  d/0 


jA;  -  r) 
A/°/  /Or  a/* 


which  is  the  same  as  equation  (10. 15)  derived  for  the  circular 
cylindrical  geometry. 


Thus,  the  solutions  to  equation  (10. 15)  also  apply  approximately  to 
equation  (10. 37)  when  the  ray  paths  do  not  approach  the  origin 
closely.  ..  - 


(1)  M.  Born  and  E.  Wolf,  Principles  of  Optics,  Pergamon,  1959, 
pp.  120—123. 

(2)  A.  L.  Nelson,  K.W.  Folley,  and  M.  Coral,  Differential  Equations, 
Heath,  1952,  pp.  44-45. 
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APPENDIX  A-ll 

RAY  PATH  GEOMETRY  FOR  BACK-SCATTERING  BY  A 
NON-UNIFORM  LOSSY  DIELECTRIC- 
COATED  METALLIC  SPHERE 

It  is  assumed  that  the  scattering  object  consists  of  a  perfectly  reflective 
(metallic)  sphere  coated  with  a  layer  of  lossy  dielectric  material.  The 
thickness  of  the  dielectric  coating  is  assumed  to  be  constant  all  around 
the  sphere.  The  electromagnetic  properties  of  the  dielectric  material 
are  assumed  to  vary  only  in  the  angular  direction  away  from  an  axis 
of  symmetry  and  not  angularly  around  the  axis  of  symmetry  nor  radially 
from  the  center  of  the  sphere.  The  ray  path  geometry  of  a  plane  electro¬ 
magnetic  wave  incident  on  these  concentric  spherical  boundaries  along 
the  axis  of  symmetry,  and  the  resultant  scattered  wave  ray  path  geo¬ 
metry,  necessary  for  the  determination  of  the  scattering  cross  section, 
can  be  completely  defined  by  the  ray  paths  in  the  plane  of  incidence 
(as  defined  in  Appendix  4).  A  normal  view  of  the  plane  of  incidence, 
showing  the  equatorial  section  of  the  spherical  boundaries  and  containing 
the  axis  of  symmetry,  is  shown  in  Figure  All-1.  The  angles  of  refraction 
and  the  curvature  of  the  ray  paths  in  the  dielectric  material  are  indica¬ 
tive  of  a  dielectric  having  a  relative  phase  factor  with  a  value  less  than 
unity  at  the  axis  of  symmetry  at  the  front  of  the  sphere  and  increasing 
toward  unity  in  the  angular  direction  away  from  this  location.  Three 
characteristic  ray  path  configurations  are  apparent,  but  only  the  con¬ 
tributions  to  the  scattered  wave  due  to  "single-transit”  paths  through 
the  dielectric  are  detailed.  It  should  be  noted  that  contributions  to  the 
scattered  wave  also  come  from  rays  that  have  been  carried  around  the 
back  of  the  sphere  in  the  dielectric  coating.  The  geometrical  details  of 
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SPHERE  SURFACE 


Figure  All-1  Ray  Path  Geometry  in  the  Plane  of  Incidence 
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each  of  the  three  characteristic  ray  paths  must  be  known  in  order  to 
determine  the  relationship  between  the  incident  wave  and  the  scattered 
wave. 

The  "  v  ”  ray  path  in  this  case  is  identical  to  the  "  v"  ”  ray  oath  for 
the  case  of  a  uniform  dielectric,  the  details  of  which  are  given  in 
Appendix  5. 

The  "  VcL  ”  ray  path,  which  must  encircle  the  sphere  at  least  once 
in  order  to  contribute  to  back- scattering,  may  be  neglected  for  thin 
coatings  (see  Appendix  5). 

Consider  the  M  Vdts”  ray  path,  in  which  the  ray  is  transmitted  through 
the  outer  dielectric  surface,  passes  through  the  dielectric  coating 
until  it  meets  the  sphere  surface,  is  reflected  from  the  sphere  surface, 
again  passes  through  the  dielectric  coating  until  it  meets  the  outer 
dielectric  surface,  and  is  transmitted  through  the  outer  dielectric 
surface  to  the  receiver  location.  In  general,  such  a  ray  path  may 
traverse  the  dielectric  coating  2  m  times  before  emerging.  This  path 
is  detailed  in  Figure  All -2.  It  is  obvious  from  this  figure  for  the 
"far-field”  condition  (  oo  )  that 

(11.1)  =  Qiyr  +  +  0,  -*■  0g, 

analogous  to  the  situation  for  a  uniform  dielectric  coating  (see 
Appendix  5).  It  is  necessary  to  express  dtv>  &\  and  in 

terms  of  Q  \,r  in  order  to  define  the  ray  path  geometry. 
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Appendix  A- 4  gives 

«i.2)  .  YT  iu  6Lv 

Ax /r\  nL  j  =  . . - . . . r. 

^6^  9i\r  +  +  j6r  |  * 

and 

«i.3)  W  9tv  =y=^^  *u  [ e* +  («?  +  ]4  ] 

where  €.y  is  the  dielectric  constant  and  is  the  relative  con¬ 
ductivity  of  the  dielectric  material. 

It  is  shown  in  Appendix  A- 10  that  when  the  sphere  radius  r  is  large 
compared  to  the  coating  thickness  A  this  spherical  geometry  can 
be  approximated  by  a  circular  cylindrical  geometry  for  small  angles 
away  from  the  axis  of  symmetry,  as  shown  in  Figure  All-3.  In  this 
coordinate  system  the  axis  of  symmetry  remains  the  same  and  the 
plasma  properties  vary  only  in  the  /°  direction  and  do  not  vary  in  the 
z.  direction  nor  angularly  around  the  axis  of  symmetry. 

The  equations  of  transformation  between  these  two  coordinate  systems 
(see  Appendix  A- 10)  are 

(11.4)  (°  s  (r +  -&)  6 
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Figure  All -3  Approximation  for  Characteristic  Ray  Path  ”vdsM 
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(11.5)  Z  s  (  r  +  A.  )  -  R. 

where  R  is  the  distance  from  the  origin  of  the  spherical  coordinate 
system  to  a  point  on  the  ray  path.  Defining  the  angle  Iff  to  be  the 
angle  between  the  positive  z-axis  (axis  of  symmetry)  and  the  tangent 
to  the  ray  path  at  any  point  (measured  in  the  circular  cylindrical 
system)  gives 

(11.6)  tom.  f  = 

oL  21 

Using  equations  (11.4)  and  (11.5)  in  equation  (11. 6)  gives 

(n.7)  w  f  =  _|v+jr,)  AL 

Defining  the  angle  X  to  be  the  angle  between  the  axis  of  symmetry 
and  the  tangent  to  the  ray  path  at  any  point  (measured  around  the  origin 
in  the  spherical  system)  gives 

(li.s)  tom.  (x-0)  =  R-fn* 

Equations  (11.7)  and  (11.8)  give 

(n.9)  W  ^  W\|-x  -fl) 
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It  can  be  seen  from  equations  (11.4),  (11.5),  and  (11.9)  that  a  ray 
path  in  the  spherical  system  meeting  the  outer  spherical  dielectric 
surface  (  R  *  r  JOc  )  at  9  *  9  i,»r  and  parallel  to  the  axis  of 
symmetry  (an  incident  wave  ray  path)  is  transformed  into  a  ray  path  in 
the  circular  cylindrical  system  meeting  the  outer  plane  dielectric 
surface  (z  =  0)  at  /°*/^  and  at  an  angle  to  the  axis  of  symmetry 

(Z-  axis)  where 

(11.10)  /%  -  |r  +  -£. )  eLir 

(11.11)  =.  0iu. 

This  is  shown  in  Figure  All -3.  It  can  also  be  seen  that  a  ray  path  in 
the  circular  cylindrical  system  leaving  the  outer  plane  dielectric 
surface  (z  =  0)  at  /°  =  /°^  and  at  an  angle  ifr^  *  IT  - 
axis  of  symmetry  (z  -  axis)  is  transformed  into  a  ray  path  in  the 
spherical  system  leaving  the  outer  spherical  dielectric  surface  (R= 
r  +&)  at  Q  *  6/^  and  at  an  angle  X^  to  the  axis  of  symmetry 
(a  scattered  wave  ray  path)  where 

(11.12) 

(ii.  is)  + 

From  equations  (11. 1),  (11.  2),  and  (11. 13)  it  can  be  seen  that  X ^  is 
actually  0  j,  R  *  Usin8  equations  (11. 3)  and  (11. 13)  gives 
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(11. 14) 


0i.R>rJ. 


<rw 


I  vx 


2.  »  »• 
6'T'  +  ^ 


i'4 


L<Lrrv 


Thus  it  is  only  necessary  to  find  and  @ieLrm.  in  terms  of 

Q  ^y.  in  order  to  define  the  scattered  wave  ray  paths. 

The  actual  ray  path  geometry  in  the  dielectric  layer  involves  m 
reflections  from  the  metallic  wall  and  m-1  reflections  from  the  inner 
surface  of  the  outer  dielectric  boundary  before  the  ray  path  emerges 
from  the  dielectric  layer  at  •  Because  there  is  no  variation  of  the 

electromagnetic  characteristics  of  the  dielectric  material  in  the  z 
direction  (parallel  to  the  axis  of  symmetry)  the  part  of  the  ray  path  from 
the  metallic  wall  back  to  the  outer  dielectric  boundary  can  effectively 
be  replaced  by  its  image  in  the  metallic  wall  in  order  to  find  /°t  and 
^LdLl  ,  as  shown  in  Figure  All -3.  This  technique  can  be  re¬ 
peated  for  all  successive  round  trips  through  the  dielectric  layer.  Thus, 
only  a  single  effective  ray  path  proceeding  in  the  positive  z  direction 
need  be  considered.  This  ray  path  starts  at  the  dielectric  surface 
(  /°»  /°o  ;  z.  =  0  )  at  some  initial  angle  to  the  axis  of  symmetry 
(  •  0tJL )  .  The  values  of  and  Q  uirm.  for  the  actual 

path  are  the  same  as  the  values  of  p  and  along  the  single 
effective  ray  path  at  the  position  defined  by  z  =  2  m  h. 
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The  electromagnetic  characteristics  of  the  dielectric  material 
(dielectric  constant,  relative  conductivity,  and  relative  phase  factor) 
vary  symmetrically  with  respect  to  6  and  consequently  each  can 
be  expressed  in  the  form  of  a  Maclaurin's  series  expansion  of  the 
form 


(u.i5)  n  =  <ia  +  e*  +  $+  &4  +  — 


which  can  also  be  written  in  the  form 


(11.16)  %  =  jf 


+  Ml 


/o*  + 


where  only  the  first  two  terms  of  the  series  need  be  considered  for 
regions  near  the  axis  of  symmetry.  In  Appendix  10  the  quadratic 
form  of  variation  of  the  relative  phase  factor  was  analyzed  and  the 
ray  paths  determined  for  the  circular  cylindrical  system.  Such  a 
ray  path,  with  the  initial  conditions  /®  =/*  ;  *  O  and 

t,  *  Oti  ,  is  the  effective  ray  path  necessary  to  determine  0(,«jlrtw 
and  (in,  from  which  ^  i.  v*  3  nrv  may  be  found  by  equation 
(11. 14).  Representing  the  functional  form  of  the  variation  of  the 
relative  phase  factor  near  the  axis  of  symmetry  by  the 

expression 


(11.17) 
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The  equation  of  the  effective  ray  path,  from  Appendix  10,  is 

[—  j  M.  [  (•frx)'/l'  zj  = 

(11.18)  * 

when  fir  is  algebraically  positive,  and 

<u-i9>  „  .  „ 
/•(“^-rt‘1  -/s  (!z!£^-'°T 

when  -&•  is  algebraically  negative,  and 
where 


(11.20) 


1  +  iom.2, 

fir*  ( I  + 


Note  that  Bti  can  be  expressed  in  terms  of  @i\r  by  means  of 
equation  (11.2).  Solving  equations  (11.18)  and  (11. 19)  for  /°  in 
terms  of  z  gives 


(11.21)  p  * <^rs-  +  K^Ai»dk  [A.,  (AK)^  z] 

+  oaA  [/3r<  (AKj^z] 
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when  A  is  algebraically  positive,  and 


(11.22)  _£  j  ^  (-AK)’4z] 

+/2  [  A,  (-4k)  1  *•] 


when  A  is  algebraically  negative.  Thus,  from  equations  (11. 21) 
and  (11. 22),  is  given  by 


(11.23) 


+  /g  CdiM.  [#•$ 


when  A  is  algebraically  positive,  and 
(n.24) 

cad.  ^  /3f$  (-AK)^*  2,/m.A.] 
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when  4-  is  algebraically  negative. 

From  equations  (11.21)  and  (11.22)  it  can  be  seen  that 


*“ ■  -  As  /©  [/As  (-4k)  z:j 

(11.25) 

+  As  (*k|a(  j’tul*  (  Ak),4z] 


when  §■  is  algebraically  positive,  and 


|-4k)  /l/«  W  [  A,  (-Ak  )  /lx ] 


(11.26) 


+  Aj  |-4k) 


Vz 


l-K~'As“2 


-4 


-/J2]l4cM.[/9I.s|-4K)!iz.] 


when  4r  is  algebraically  negative.  Thus,  from  equations  (11.6), 
(11.25),  (11.26),  and  the  previously  stated  fact  that  Is 

equal  to  ,  tan  is  given  by 
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(11.  27) 

Wt  Siirrn_= /3f}  (lk)/2/f  Avn-ft. [#r,  Uk)  ^2  <*«•&■ 

+a5  (4k)\^ 

when  &  is  algebraically  positive,  and 
(11.28) 

la'".  "A,  l-U)^  Ain.  [a,  (-•fcl<)^'2nn.ft. 

+  A,  [/<?rs  (-AKj^S/mil] 

when  4r  is  algebraically  negative.  Thus,  equations  (11. 14),  (11.23), 
(11.24),  (11.27),  and  (11.28)  give  ^i,R>rd.Smv  as 
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(n.29)  jTTjj-  c*^[A«(^r2/m-*-J 


,“i  .  -2. 


+  axe  M n  (  -L  (An  o^tSUv 

VF 


Ai*-4.  j  firt  (4K|  /Z2 


/®r.  l*<lVa 


+  /3rs  |<K)4[  ♦/o«j4c*4  !  ^ 


1 


1 

+  €**  +  ffj1  CeA  ^  axe 'Em* 


r 


_  , .  ,‘/z 

^r,l<K|  /°o 


/* 


**4.1  /8„UK|  u«. 


+^s  1**1 


*4  /  -i  —  -x  !4 

^,‘k 


AjlK)  Zmi. 


rs 


i 

j  L 


when  *w  is  algebraically  positive,  and 
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(11.30)  9, 


Vt 


IR  vrds  on 


Ctri 


2/vn.&j 


f*(r  +  *)*  (v**| 


2.  y 


+  dRC^AV 


1 

«  . . r  ( 

/ 

IT 

>4uv,  OJl£.ta/n. 

vr 

> 

-^,(-WyVe 


AWt 


rs 


l  -+ 


-i  „  -* 


!-K  /S„ 

o 


'A 


”  -/*«  i 


/^.s("4k)  ZmR. 


£  -  a  -a 


€r  +  I  <Er  +■  (T^.  Cod.  QActayvw 

^k,(-4K)/£Zmv^J 

& 


-^Hk)  /> 


+  /,WW’liZ^^om^ 


:a 


/3ri(-Akr)  2mvJt 


•4 


when  4  is  algebraically  negative. 
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Thus,  for  ray  paths  near  the  axis  of  symmetry,  and  an  incident 
wave  propagating  in  a  direction  parallel  to  the  axis  of  symmetry, 
the  ray  path  geometry  is  defined  for  an  axially- symmetric  non- 
uniform  lossy  dielectric  coating  on  a  metallic  sphere  when 
/O  /(\r-*-&)  is  replaced  by  9 ir  and  ir  is  replaced  by 
2A(e/lr+A)  in  the  preceding  equations. 
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APPENDIX  A-12 

DIVERGENCE  FACTOR  FOR  BACK- SCATTERING  BY  A 
NON-UNIFORM  LOSSY  DIELECTRIC-COATED  METALLIC 

SPHERE 

It  is  assumed  that  the  scattering  object  consists  of  a  perfectly  reflective 
(metallic)  sphere  coated  with  a  layer  of  lossy  dielectric  material.  The 
thickness  of  the  dielectric  coating  is  assumed  to  be  constant  all  around 
the  sphere.  The  electromagnetic  properties  of  the  dielectric  material 
are  assumed  to  vary  only  in  the  angular  direction  away  from  an  axis 
of  symmetry  and  not  angularly  around  the  axis  of  symmetry  nor 
radially  from  the  center  of  the  sphere. 

The  geometrical  divergence  factor  for  a  plane  electromagnetic  wave 
propagating  in  a  direction  parallel  to  the  axis  of  symmetry  and  inci¬ 
dent  on  the  non-uniform  lossy  dielectric-coated  metallic  sphere, 
necessary  for  the  determination  of  the  scattering  cross  section,  is 
determined  by  the  ray  path  geometry  (as  derived  in  Appendix  A- 11). 

The  general  form  of  the  geometrical  divergence  factor  for  back- 
scattering  from  concentric  spherical  regions,  from  Appendix  A-6,  is 


(12. 1)  DiR 


2  Rr  jSOm.  AfljR 

A0i(r-o  */**■  A  6ir 


where  it  is  necessary  to  express  &  iR  explicitly  in  terms  of 
6 1 r  before  the  limit  can  be  taken*  Equation  (12. 1)  can  also  be 
written  in  the  form 


A12-1 


TR63-217A 


(12.2) 


4Ri 


J Qionrv 


°iR  =  |t+a)z  A$ir-~o  rzzrze 17 


The  two  different  functional  relationships  between  and  &i\r 

are  given  in  Appendix  A- 11  for  the  two  different  characteristics  ray 
paths  that  contribute  to  the  back- scattering  cross  section,  thus 
leading  to  different  divergence  factors  for  the  two  different 
characteristic  ray  paths. 


The  back- scattering  geometrical  divergence  factor  for  the  ”  vr  ” 
ray  path  in  this  case  is  the  same  as  for  the  "  \r  ”  ray  path  in  the 
case  of  a  uniform  dielectric,  and  is  given  ini  Appendix  A-6. 

The  back- scattering  geometrical  divergence  factor  for  the  "  vis  " 
ray  path  is  given  by  two  different  equations  for  the  two  different  forms 
of  variation  of  the  dielectric  properties  near  the  axis  of  symmetry 
(see  Appendix  A-ll).  Substituting  the  formulas  for  0  iR  \rJLStm. 
from  Appendix  A-ll  into  equation  (12. 2)  and  taking  the  limit  gives 


(12.3) 


4Rr 


aVZ~  <m.A. 


I 

+  z  vr 


KVt  h2-/9-*  */s] 


\ 


'/z 


2  VST  mv 

TTTI 


•Mr 

nil 
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when  ^y0  is  algebraically  positive,  and 

4R 


(12.4)  = 


( 


l  r r. .  i,''*  aVF mu] 

11*  l  l>  |r.«|  J 


Z-yfz"  (""M^  [l+z  A*  -^/a] 


^  [|-V|/Z 


2.VjT  onfl. 
|r+4l) 


r 


when  V  is  algebraically  negative. 
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APPENDIX  A- 13 

GEOMETRICAL  OPTICS  ELECTROMAGNETIC  WAVE 
BACK-SCATTERING  CROSS  SECTION  OF  A  NON-UNIFORM 
LOSSY  DIELECTRIC-COATED  METALL I C  SPHERE 

Detailed  here  is  the  theoretical  determination  of  the  electromagnetic 
wave  back- scattering  cross  section  (radar  cross  section)  of  a  per¬ 
fectly  reflective  (metallic)  sphere  coated  with  a  non-uniform  layer 
of  lossy  dielectric  material  via  the  geometrical  optics  approach  using 
the  ray  path  field  intensity  addition  technique  (see  Appendix  A-3).  The 
thickness  of  the  dielectric  coating  is  assumed  to  be  constant  all 
around  the  sphere.  The  electromagnetic  properties  of  the  dielectric 
material  are  assumed  to  vary  only  in  the  angular  direction  away 
from  an  axis  of  symmetry  and  not  angularly  around  the  axis  of 
symmetry  nor  radially  from  the  center  of  the  sphere.  The  incident 
electromagnetic  wave  is  taken  to  be  propagating  in  a  direction  parallel* 
to  the  axis  of  symmetry. 

The  only  difference  between  the  present  situation  with  a  non-uniform 
dielectric  and  the  situation  with  a  uniform  dielectric  is  the  ray  path 
geometry,  and  consequently  the  geometrical  divergence  factors  for 
the  various  ray  paths  involved  are  different.  The  results  for  the 
uniform  dielectric  situation  given  in  Appendix  A-7  can  be  used  for  the 
non-uniform  dielectric  situation  simply  by  substitution  of  the 
appropriate  divergence  factors  into  the  equations.  Thus,  using  the 
results  from  Appendixes  A-7  and  A- 12,  the  normalized  radar  cross 
section  is  given  by 
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(13.1) 


nr 


=  i 


*1 


III *r£) 

fsr* i 


I 


OTi- I 


P\r&  ^  € 


where 


(13.2) 

hmn 


o^M.  U/2 


zVz~  omJL 

l»+*-)  . 


+ 


zViT 


[  1  ■*■  z  /®C,  ^/«] 


v*  -v  A 


-I 


when  is  algebraically  positive,  and 


(13.3) 


Dm,- 


+  zVT"  I"  M  [i  +  2/S^  A*] 


(— -%a) 


'4  zVz~  mvd 


|r**)  J 


-I 
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when  k/6  is  algebraically  negative.  It  can  be  seen  that  when 
•ft  £  is  zero  (uniform  dielectric)  equations  (13.  2)  and  (13. 3)  both 

reduce  to  the  form 


.  ft  I  I"1 

1  +  ^  [77*1  atJ 

This  is  not  exactly  the  same  as  the  correct  form  of  as 

originally  derived  for  the  uniform  dielectric  situation,  because  the 
ray  path  equation  for  the  spherical  geometry  was  approximated  in 
order  to  solve  it  for  the  non-uniform  dielectric  situation.  In  order  to 
preserve  the  correct  asympotic  form  of  it  is  necessary  to 

replace  (r  +  h)  by  r,  which  is  essentially  one  of  the  approximations 
that  was  made  in  simplifying  the  spherical  geometry  problem.  Thus, 
the  equation  for  the  relative  geometrical  divergence  factor  becomes 


(13.4) 


'rm.  - 


(13.5) 


(V)14 


+  ("K« )  [l  +  Z/3^  ^a] 


IVrJk  \  (ft/»)!4  ^ 
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when  V  is  algebraically  positive,  and 

!4  aVzT  iw.fl  1 

+  zTe”  /^l‘s  ['  z'9**  V] 

^  [  | -^)*  «y*  ^1 

when  ^  is  algebraically  negative .  Note  that  the  values  of  the 
dielectric  properties  to  be  used  in  the  above  equations  are  the 
values  at  the  axis  of  symmetry. 

Some  insight  into  the  behavior  of  the  radar  cross  section  as  given  by 
equations  (13. 1),  (13.5),  and  (13.6)  can  be  obtained  by  considering 
the  individual  terms  that  contribute  to  the  overall  result.  Consider 
a  thin  pencil  beam  of  electromagnetic  radiation  incident  along  the 
axis  of  symmetry  on  the  non-uniform  lossy  dielectric-coated  metallic 
sphere.  At  the  outer  dielectric  surface  part  of  the  incident  power  is 
reflected  in  a  diverging  cone.  This  cone  is  the  same  for  any  spatial 
variation  of  the  dielectric  properties  within  the  dielectric  layer,  and 
is  shown  in  Figure  A13-1  for  a  uniform  dielectric,  a  dielectric  with 
a  de-focusing  gradient  (  positive),  and  a  dielectric  with  a  focusing 
gradient  ( negative).  This  contribution  to  the  scattered  wave 
power  is  represented  by  ,  the  first  term  in  the  brackets  of 

equation  (13. 1).  Subsequent  terms  (the  infinite  sum)  in  equation 
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(13. 1)  represent  the  contributions  from  the  part  of  the  incident 
power  that  is  transmitted  into  the  dielectric  layer  and  subsequently 
is  partly  transmitted  out  again  after  each  successive  round  trip 
through  the  dielectric  layer.  The  emergent  cones  after  the  first 
and  second  round  trips  through  the  dielectric  coating  are  shown  in 
Figure  A13-1,  for  the  three  different  forms  of  spatial  variation  of 
the  dielectric  properties  within  the  dielectric  layer.  It  can  be  seen 
that  when  a  de-focusing  gradient  is  present  the  emergent  cone  is 
wider  and  consequently  the  power  density  is  smaller,  and  vice- 
versa  when  a  focusing  gradient  is  present.  It  is  important  to  remem¬ 
ber  that  the  relative  phasing  between  the  electric  fields  associated 
with  each  of  the  scattered  cones  of  power  determines  the  resultant 
total  scattered  power  density.  Hence,  relative  to  the  uniform  dielectric 
layer,  the  dielectric  layers  with  either  focusing  or  de-focusing 
gradients  can  produce  either  increased  or  decreased  total  back- 
scattered  power  density.  However,  for  very  thin  dielectric  layers  the 
usual  effect  of  a  de-focusing  gradient  is  a  relative  decrease  in  the 
radar  cross  section  and  vice-versa  for  the  focusing  gradient.  As  in 
the  case  pf  a  uniform  dielectric  coating,  conditions  can  be  found  for 
which  the  radar  cross  section  vanishes.  Also,  equation  (13. 6)  shows 
that  when 

(13.7)  +  -1=  fij 
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the  radar  cross  section  becomes  infinitely  large.  This  is  because  the 

1L 

rn  emergent  cone  of  power  is  focused  to  a  point  at  the  radar  receiver 
location,  producing  an  infinitely  large  power  density  there.  This  result 
is  not  valid  because  one  of  the  applicability  criteria  for  the  geomet¬ 
rical  optics  solution  has  been  violated  (the  solution  is  not  good  in  the 
vicinity  of  a  focal  point).  However,  this  condition  is  indicative  of  the 
fact  that  the  back- scattered  power  density  will  be  very  large  and 
consequently  significant  increases  in  the  radar  cross  section  may  be 
observed  under  these  conditions. 
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APPENDIX  A-14 

SPATIAL  VARIATION  OF  PLASMA  PROPERTIES  ' 

IN  THE  SHOCK  LAYER  AROUND  THE  NOSE  OF  A  SPHERE  IN 
HYPERSONIC  FLIGHT 

The  computation  of  the  electromagnetic  properties  of  the  ionized  gas 
(plasma)  in  the  shock  layer  around  the  nose  of  a  blunt  body  such  as  a 
sphere  in  hypersonic  flight  thru  air  has  received  considerable  atten¬ 
tion  in  recent  years.  The  problem  must  be  handled  numerically  and 
is,  in  general,  very  complicated. 

When  the  aerodynamic  parameters  are  such  that  chemical  equilibrium 
is  maintained  throughout  the  flow  field  in  the  nose  region, the  situation 
is  greatly  simplified.  The  dependence  of  the  plasma  properties  on  the 
absolute  size  of  the  body  practically  disappears  and  consequently 
normalized  solutions  can  be  obtained.  The  aerodynamic  variables 
involved  in  this  situation  are  the  composition,  temperature,  and  density 
of  the  ambient  air  thru  which  the  body  flies,  and  the  speed  of  the  body. 

-  Charts  of  the  flow  field  properties  of  interest  for  electromagnetic 
considerations  have  been  prepared^  for  special  locations  in  the 
flow  field,  such  as  the  stagnation  point.  A  typical  example  of  the 
dependence  of  the  plasma  properties  at  the  stagnation  point  on  the 
speed  of  the  body  is  shown  in  Figure  A14-1,  derived  from  Reference  1. 
Both  the  plasma  frequency  and  the  electron  collision  frequency  are 
shown,  normalized  to  a  reference  frequency  of  35Gc  for  an  ambient 
air  pressure  of  10  mm  Hg,  which  is  approximately  equivalent  to  an 
altitude  of  100,000  feet.  The  spatial  variation  of  the  plasma  proper¬ 
ties  in  the  shock  layer  around  the  body  away  from  the  stagnation 
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Figure  A14-1  Stagnation  Point  Plasma  Properties  for  Equilibrium 
Conditions 
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point  depends  on  the  initial  conditions  at  the  stagnation  point  and  the 
exact  shape  of  the  blunt  nose  of  the  body. 

For  a  spherical  nose,  the  pressure  variation  in  the  shock  layer  from 
the  stagnation  point  to  45  degrees  off  axis  approximately  follows  the 
Newtonian  formula 

=  cetJ-fr 


(14. 1) 


where  p  is  the  pressure  in  the  shock  layer,  p  is  the  pressure  at  the 
stagnation  point,  and  Q  is  the  angle  measured  around  the  sphere 
from  the  stagnation  point.  The  electron  collision  frequency  in  air 
is  almost  directly  proportional  to  the  pressure,  ’  hence  the 
variation  of  the  electron  collision  frequency  \)  around  the  nose  of 
the  sphere  is  given  by 


where  ~\)  is  the  value  of  the  electron  collision  frequency  at  the 
stagnation  point.  The  expansion  of  the  hot  compressed  air  from  the 
stagnation  region  as  it  flows  around  the  body  approximately  follows 
an  isentropic  process.  Reference  to  charts  of  high-temperature  air 
properties  ^  shows  that  for  isentropic  expansions  over 

small  pressure  changes  in  the  parameter  regime  shown  in  Figure 
A1 4-1,  the  plasma  frequency  decreases  approximately 

proportionally  to  the  pressure,  hence 
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(14.3) 


jfe-  £  coi,1 6 


where  f  is  the  value  of  the  plasma  frequency  at  the  stagnation 
ps 

point.  Near  the  stagnation  point,  equations  (14.  2)  and  (14. 3)  can  be 
written  in  the  forms 


(14.4)  2  t-fi*1 

and 


(14.5) 

=  »-  0** 

In  general,  for  a  less  restrictive  situation,  the  plasma  frequency 

and  the  electron  collision  frequency  in  the  vicinity  of  the  stagnation 

point  can  be  expressed  in  the  forms 

(14.6) 

V 

£  i  ♦  fe*.  6* 

and 

(14. 1) 

i 

.a  \  + KB* 
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Using  the  normalized  plasma  frequency  (  XX  -  "ff©  / *f  )  and 
the  normalized  electron  collision  frequency  (  XI  c  *  ) 

these  expressions  become 

(14.8)  -fa.  =  XI  ^ 
and 

ne  =  ncs  |i  +  *c  0*1 

(14.9) 


where  is  the  angular  frequency  of  the  incident  electromagnetic 
wave.  The  electromagnetic  wave  relative  phase  factor  for 
plasma,  introduced  in  Appendix  1,  is  given  by 


(14.10)  ^  = 
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Substituting  equations  (14. 8)  and  (14. 9)  into  (14. 10)  and  neglecting 
powers  of  9  greater  than  the  second  gives 


(14.  ll)  /3r  S  0*) 


where 


(14.12) 


and 
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(14.13) 

v  ^ 


l+a<is 


The  variation  of  the  relative  phase  factor  of  the  plasma  layer  near 
the  stagnation  point  is  important  in  the  determination  of  the  nose-on 
radar  cross  section  of  the  plasma  covered  blunt  body.  Both  the 
plasma  frequency  and  the  electron  collision  frequency  decrease  in 
going  away  from  the  stagnation  point  in  flow  fields  around  blunt 
bodies  (fe^and  are  algebraically  negative).  This  is  depicted  in 
Figure  A14-2.  Equations  (14. 11)  thru  (14. 13)  show,  however,  that 
very  near  the  stagnation  point  the  relative  phase  factor  may  either 
increase  or  decrease  at  first  in  going  away  from  the  stagnation  point 
(  kp  may  be  either  positive  or  negative  algebraically).  This  is  also 
depicted,  in  exaggerated  form,  in  Figure  A14-2.  For  example,  when 
there  is  no  angular  variation  of  the  electron  collision  frequency 
(  k  t  *  O  )  the  variation  of  the  relative  phase  factor  is  opposite 
to  that  of  the  plasma  frequency  when  the  normalized  plasma  frequency 
is  less  than  VzT  ,  and  in  the  same  direction  when  the  normalized 
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Figure  A14-2  Typical  Variation  of  Equilibrium  Plasma  Properties 
Around  a  Sphere 
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plasma  frequency  is  greater  than  VzT  ,  that  is 

(14. 14)  i-k^i  =  -***.(4^)  {^;<0vr 

(14.15)  L~>V5r 

When  the  relative  phase  factor  increases  monotonic  ally  in  going  away 
from  the  stagnation  point  an  incident  electromagnetic  wave  is  de- 
iocused  relative  to  the  back-scattering  direction,  as  discussed  in 
Appendix  13.  When  the  relative  phase  factor  decreases  at  first  in 
going  away  from  the  stagnation  point  the  part  of  an  incident  electro¬ 
magnetic  wave  within  the  region  of  decreasing  relative  phase  factor 
is  focused  in  the  back-scattering  direction.  Thus,  in  general,  de- 
focusing  occurs  at  low  values  of  plasma  frequency  and  focusing 
occurs  at  high  values  of  plasma  frequency. 

The  prediction  of  the  spatial  variation  of  the  plasma  properties  in  the 
flow  field  when  chemical  equilibrium  is  not  maintained  throughout  the 
flow  field  is  extremely  difficult.  The  composition  of  the  gas  through¬ 
out  the  flow  field  will  depend  on  the  absolute  size  of  the  body  as  well 
as  its  generic  shape  because  the  chemical  reactions  proceed  at  finite 
rates  and  the  flow  time  around  the  body  is  proportional  to  the  body 
size.  Thus,  a  very  large  body  gives  long  flow  times  and  consequently 
near -equilibrium  conditions  can  be  maintained  throughout  the 
shock  layer  around  the  body.  A  very  small  body  gives  short  flow 
times,  the  flow  field  chemical  reactions  will  be  far  from  equilibrium, 
and  consequently  the  flow  field  properties  will  depart  significantly 
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(5) 

from  the  equivalent  equilibrium  properties.  A  typical  calculation'  ’  of 
the  non- equilibrium  flow  field  properties  around  the  nose  of  a 
sphere  in  hypersonic  flight  thru  air  is  shown  in  Figure  A14-3.  It  can 
be  seen  from  this  figure  that  the  electron  collision  frequency  does 
not  change  very  much  thru  the  shock  layer  from  near  the  shock  front 
to  near  the  body  (in  the  radial  direction),  and  has  a  slight  decrease 
in  the  angular  direction  away  from  the  stagnation  point.  Actually, 
within  the  angular  range  shown  (from  zero  to  approximately  20 
degrees)  the  variation  of  the  electron  collision  frequency  follows  the 
form 


(14.16)  -O-C.  -  ftcs 

which  is  the  same  as  for  the  equilibrium  flow  field.  The  plasma  fre¬ 
quency  shows  significant  variation  in  both  the  radial  direction  thru  the 
shock  layer  and  the  angular  direction  away  from  the  stagnation  point. 
Within  the  angular  range  shown,  the  variation  of  the  plasma  frequency 
in  the  angular  direction  follows  the  form 

U4.17)  1 1  ♦  e *■ ) 


where  W^is  between  -3  and  -5,  which  gives  a  much  more  rapid  change 
than  predicted  for  the  equilibrium  flow  field. 

The  shock  detachment  distance  at  the  stagnation  point,  for  the  same 
conditions  as  given  in  Figure  A14-1,  is  shown  in  Figure  A14-4.  As 
usual,  this  shock  layer  thickness  (h)  is  small  compared  to  the  nose 
radius  (r)  of  the  body. 
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Figure  A14-4  Shock  Detachment  Distance  at  Stagnation  Point 
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APPENDIX 

GEOMETRICAL  OPTICS  NOSE-ON  RADAR  CROSS 
SECTION  OF  A  METALLIC  SPHERE  COATED  WITH  AN 
ANGULARLY  NON-UNIFORM  LOSSY  PLASMA 


Given  here  is  the  radar  cross  section  d*  ,  obtained  by  the  geometrical 
optics  approach,  of  a  perfectly  reflective  (metallic)  sphere  of  radius 
r  covered  with  a  non-uniform  coating  of  lossy  plasma.  The  plasma  is 
assumed  to  have  uniform  thickness  -K  all  around  the  sphere.  The 
electromagnetic  properties  of  the  plasma  are  assumed  to  vary  only  in 
the  angular  direction  around  the  sphere  away  from  an  axis  of  symmetry 
and  not  angularly  around  the  axis  of  symmetry  nor  radially  from  the 
center  of  the  sphere.  The  incident  electromagnetic  wave  is  taken  to  be 
propagating  in  a  direction  parallel  to  the  axis  of  symmetry  (incident 
nose-on).  Furthermore,  it  is  assumed  that  the  plasma  can  be 
characterized  (electromagnetic ally)  by  a  permeability  ’a 

dielectric  constant  ,  and  a  conductivity  <S^  analogous  to  the 

situation  for  non-ionized  media.  Thus,  the  situation  becomes  that  of 
a  non-uniform  lossy  dielectric-coated  metallic  sphere.  Using  the 
results  from  Appendixes  1,  13,  and  14,  the  normalized  nose-on 
radar  cross  section  is  given  by 


where  is  the  reflection  coefficient  for  a  plane  transverse 

electromagnetic  wave  propagating  normally  from  free  space  into  a 
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plane  interface  between  free  space  and  a  uniform  plasma  having  the 
electromagnetic  properties  of  the  coating  at  the  axis  of  symmetry, 
"tps  is  the  propagation  factor  for  a  plane  transverse  electromagnetic 
wave  propagating  thru  an  unbounded  uniform  plasma  having  the 
electromagnetic  properties  of  the  coating  at  the  axis  of  symmetry, 
and  Dm  is  the  relative  geometrical  divergence  factor  for  the  cmr  ray 
path  thru  the  plasma  coating.  Appendixes  1,  4,  13  and  14  give 


1-  t- 


(15.2)  r*s  = 


i- jn 


\x/z 


cs 


1*1- 


i-jn 


cs 


(i5.3)  = 


j4JL Vt 


i-jn 


cs 


= 


(15.4) 


Z'y/z 


cmu  [  v4  ] 

<  k'/'  [i  -  ***[*£ 


-I 


when  is  algebraically  positive,  and 


(15.5) 
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when 


is  algebraically  negative,  where 


and  in  the  vicinity  of  the  axis  of  symmetry 


(15.7) 


* + 


)+  J 

'A. 

-  - 

*.  ' 

when  the  normalized  plasma  frequency  and  the  normalized 

electron  collision  frequency  £LC  in  the  vicinity  of  the  axis  of 
symmetry  are  given  by 

us.#)  s  af9  ( •  ♦  -fcj.  6 l) 
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(i5.9)  at  a  ncs  (i  +  e‘) 


where 

I2ia  =  0d±/  co  =  normalized  plasma  frequency  (subscript  "s" 
~  indicates  "at  axis  of  symmetry") 


nc  =  /  to 

= 

co  = 


=  normalized  electron  collision  frequency 

(subscript  "s"indicates  "at  axis  of  symmetry") 

angular  plasma  frequency  (  ■2.TT  -f  ^  ) 

electron  collision  frequency 

angular  radar  frequency  (  <LTT  -f  ) 


X 


\r 


free-space  wavelength  of  the  radar  wave 


e 


angle  me  asured  around  sphere  away  from  axis 
of  symmetry. 


It  can  be  shown,  as  in  the  case  of  a  uniform  plasma  coating,  that  when 
the  plasma  coating  is  highly  reflective  1 1  r+»  I  3  )  .. 

the  radar  cross  section  <S  approaches  the  value  TTrx  (l  +  ^/r)  » 
which  is  the  value  for  a  perfectly  reflective  sphere  having  a  radius  equal 
to  the  outer  radius  of  the  plasma  coating.  When  the  plasma  coating  is 
negligibly  reflective  (  j  1"^  |  =  0  )  the  radar  cross  section 

approaches  the  value  TT  f  *-  ,  which  is  the  value  for  the  bare 

metallic  sphere,  as  would  be  expected.  The  radar  cross  section  can 
also  exhibit  maxima  and  minima  above  and  below  these  values,  as 
discussed  in  Appendix  A-13.  Specifically,  when 
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the  radar  cross  section  will  become  vanishingly  small,  and  when 


the  radar  cross  section  will  become  infinitely  large  in  the  geometrical 
optics  approximation.  This  latter  condition  can  occur  only  when 
is  algebraically  negative,  which  is  the  focusing  regime.  As  pointed 
out  earlier,  the  numerical  results  of  the  geometrical  optics  approach 
are  not  valid  under  this  condition. 
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theoretical  models  of  electromagnetic  wave  theoretical  models  of  electromagnetic  wave 

scattering  from  plasma-covered  bodies  are  scattering  from  plasma-covered  bodies  are 

formulated  and  analyzed.  The  theoretical  formulated  and  analyzed.  The  theoretical 

results  of  these  analyses  show  both  absorption  results  of  these  analyses  show  both  absorption 
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